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1.  GENERAL  OBJECTIVES 


The  research  under  contract  DAJA45-93-C-0046  and  its 
modifications  was  aimed  at  the  study  of  several  mathematical 
models  for  ground  freezing  processes. 


The  research  started  in  September  1993  with  Prof.  Mario 
Primicerio  as  the  principal  investigator.  During  1995  Prof. 
Primicerio  had  to  quit  the  team,  following  his  election  to 
Mayor  of  Florence  and  he  was  replaced  by  Prof.  Antonio  Fasano. 


The  research  has  been 
Dr •  Y .  Nakano  at  CRREL  and  i 
Federico  Talamucci.  The  last 
March/April  1996  and  was 
conclusions  and  in  pointing 
further  investigation. 


carried  out  in  close  contact  with 
ncluded  two  visits  at  CRREL  by  Dr. 
one  took  place  during  one  week  in 
very  useful  in  drawing  some 
out  further  subjects  which  deserve 


We  can  summarize  the  specific  objects  of  the  research 
under  the  following  statements: 

(a)  reviewing  the  extant  models  on  ground  freezing  according  to 
the  classification  in  (i)  models  with  a  sharp  interface  between 
frozen  and  unfrozen  region,  (i i)  models  with  a  frozen  fringe 

(b)  examining  the  possibility  of  describing  heaving  with  or 
Without  ice  lens  formation  by  means  of  models  of  class  (i) 

(c)  studying  models  of  class  (ii)  under  various  degrees  of 
approximation . 


Models  of  type 
frost  heave  models” 9 
(class  (ii))  will  be 


(i)  will  be  also  referred  to  as 
while  ground  freezing  with  frozen 
called  “secondary  frost  heave”. 


“primary 

fringes 


In  the  first  stage  it  also  included  a  study  „„  thawing 
[l].  The  work  on  classification  of  models  has  been  very  large 
and  is  incorporated  in  Talamueci's  Doctoral  Thesis  [2],  It  will 

be  used  to  produce  a  survey  paper  to  be  submitted  to  Surveys  of 
Mathematics  for  Industry. 


2.  PRIMARY  FROST  HEAVE  MODELS 


The  study  of  models  with  a  sharp  interface  between  frozen 
and  unfrozen  region  has  been  undertaken  under  an  innovative 
point  of  view.  It  has  to  be  stressed  that  not  all  the  models  of 
this  kind  presented  in  the  literature  can  be  considered 
thermodynamically  consistent. 

In  [3]  a  correct  model  has  been  selected  and  treated  under 
the  assumption  of  quasi-steady  temperature  field,  showing  that 
it  can  account  both  for  frost  penetration  and  for  lensing.  The 
complete  evolution  problem  with  ice  lens  formation  has  been 
studied  in  [4],  while  the  frost  penetration  case  has  been 
considered  in  [5]  .  In  [6]  a  numerical  work  has  been  developed 
dealing  with  the  situation  described  in  [3]  and  showing  the 
possible  occurrence  of  frost  penetration  and  of  ice  lensing. 


3.  SECONDARY  FROST  HEAVE  MODELS 


Then  the  research  has  been  oriented  towards  the  models 
with  a  frozen  fringe,  including  the  purely  thermal - i nduced 
water  flow  mechanism.  whose  existence  was  experimentally 
pointed  out  at  CRREL  by  Y.  Nakano . 

The  amount  of  work  done  in  this  direction  is  quite 
remarkable,  in  particular  considering  the  extreme  complexity  of 
the  mathematical  structure  of  the  problem.  We  quote  the  study 
of  quasi-steady  solutions  [7],  [8],  and  a  massive  series  of 
results  concerning  the  cases  of  constant  and  time  varying 
boundary  date.  It  is  particularly  remarkable  that  it  has  been 
possible  to  give  criteria  to  predict  whether  ice  lensing  or 
frost  penetration  is  going  to  take  place. 

Part  of  these  results  have  been 


comuni cated  by  Y.  Nakano 


and  F.  Talamucci  at  the  EUROMECH  Conference  333,  held  in 
Montecatini,  near  Florence,  in  June  1995  and  chaired  by  A. 
Fasano . 

The  results  achieved  not  only  provide  a  sound  mathematical 
basis  to  the  frost  heave  model  in  which  the  water  flow  is 
driven  by  the  complex  thermodynamical  evolution  of  the  frozen 
fringe,  but  are  also  a  useful  tool  to  investigate  the 
qualitative  properties  of  the  phenomenon. 

The  enclosed  report  by  F.  Talamucci  entitled  "A  quasi¬ 
steady  model  for  secondary  frost  heave  in  freezing  soils” 
contains  a  detailed  description  of  the  frost  heave  model  with  a 
frozen  fringe,  encompassing  the  fundamental  law  stating  that 
the  flux  intensity  of  migrating  water  is  the  sum  of  a  Darcy- 
like  term  (i.e.  the  usual  pressure  driven  flow)  and  of  a  term 
whose  contribution  is  a  flow  in  the  direction  opposite  to  the 
“thermal  gradient:. 


This  very  large  work  (almost  100  pages)  is  organized  in 
four  sections.  Sections  1  and  2  resume  the  main  results  already 
communicated  in  the  previous  interim  reports  and  presented  at 
the  EUROMECH  Conference,  namely:  the  complete  description  of 
the  mathematical  model  and  the  study  of  the  case  in  which  a 
constant  heat  flux  is  prescribed  at  the  boundaries  of  the 

porous  layer.  Sections  3  and  4  present  some  new  material,  i.e. 
the  analysis  of  the  time  dependent  data,  prescribed  either  as 
heat  -flux  or  as  boundary  temperatures. 

In  all  cases  existence  and  uniqueness  of  the  corresponding 
problem  is  proved  and  particular  emphasis  is  given  to  the 

possible  occurrence  of  the  ice  lensing  phenomenon.  From  the 
mathematical  point  of  view  we  have  a  set  of  partial 

differential  equations  for  pressure,  temperature,  water  flux  in 
each  of  the  three  phases:  the  frozen  and  the  unfrozen  region 
and  the  intermediate  layer  in  which  ice  and  water  coexist,  i.e. 
the  frozen  fringe.  The  interfaces  between  these  regions  are 

free  boundaries  (i.e.  they  are  unknown)  and  the  peculiar 


behaviour  of  the 

problem  is 

main  1 ; 

y  due 

to  the  fact  that 

the 

quantities  listed 

above  have 

to  satisfy 

some  constraints. 

The 

constraints  are 

un i 1 ateral 

and 

they 

are  responsible 

for 

discriminating  between  frost  penetration  and  ice  lens  growth. 

All  this  material  will  be  reorganized  in  a  series  of 
scientific  papers  for  publication  in  the  next  future. 


4.  CONCLUSIONS  AND  FURTHER  PERSPECTIVES 


As  a  conclusion,  we  can  say  that  the  work  performed  under 
the  contract  has  produced  a  significant  step  forward  in 
understanding  and  modelling  the  frost  heave  phenomenon,  taking 
into  account  some  new  factors  like  e.g.  the  thermal  driven 
water  flux  in  the  frozen  fringe  and  the  possible  rearrangements 
of  the  grains  of  the  porous  matrix  during  freezing,  generating 
a  change  of  porosity. 


From  the  mathematical  point  of  view,  a  series  of  very 
difficult  free  boundary  problems  has  been  solved,  providing  in 
some  cases  also  numerical  results.  Particularly  valuable  is  the 
qualitative  analysis,  leading  to  the  prediction  of  the  type  of 
evolution  of  the  freezing  process. 


Among  the 
cooperation  can 
distribution  of 
research  area 
material  both  f 
of  view. 


new  problems  to  be  investigated,  if  some 
be  continued,  is  a  finer  analysis  of  the 
thermal  energy  within  the  frozen  fringe.  This 
can  be  still  a  source  of  very  interesting 
rom  the  experimental  and  the  theoretical  point 
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Abstract.  The  subject  of  this  paper  is  the  mathematical  investigation  of  a  one-dimensional  secondary 
frost  heave  model  for  a  saturated  soil  with  a  frozen  fringe.  The  water  flux  qw  is  driven  by  both  thermal 
and  hydraulic  gradients;  the  volumetric  ice  content  v  is  a  known  decreasing  function  of  the 
temperature;  the  ice  in  the  frozen  fringe  is  at  rest  and  at  the  atmospherice  pressure;  the  occurring  of  the 
ice  segregation  process  is  determined  by  the  value  of  the  water  pressure  in  the  frozen  fringe. 

The  mathematical  analysis  of  the  model  has  been  developed  assuming  that  the  temperature  profile  is 
linear  in  each  region  of  the  soil  (quasi-steady  approach)  and  that  the  porous  matrix  of  the  frozen  fringe 
is  undeformable.  Besides  the  formation  of  a  pure  ice  layer,  we  consider  the  possibility  that  the  freezing 
front  moves  through  the  soil  without  a  macroscopic  break  of  the  porous  matrix  (frost  penetration). 

In  sect.  1  and  2  we  outline  and  discuss  the  equations  of  the  model  and  we  examine  the  case  in  which 
constant  thermal  fluxes  Oq,  at  the  base  and  on  the  top  of  the  soil  are  prescribed.  The  main  result  is 
concerned  with  the  possibility  of  predicting  the  kind  of  process  (lens  formation,  frost  penetration  or 
melting)  which  will  take  place  on  the  basis  of  the  knowledge  of  the  data  a0,  and  of  the  properties  of 
the  soil  (thermal  conductivities,  densities,  functions  KVK2,  v).  A  suitable  way  to  present  the  result  is 
to  display  in  a  (a^a^-plane  the  regions  where  lens  formation  or  frost  penetration  occur. 

In  sect.  3  we  consider  the  general  case  when  the  assigned  thermal  fluxes  a0,  depend  on  time.  The 
analysis  of  the  model  leads  us  to  conclude  that  the  change  in  time  of  the  boundary  fluxes  may  give  rise 
to  the  formation  of  several  lenses,  in  agreement  with  experimental  observations.  Indeed,  if  a0,  a1  verify 
certain  conditions  involving  the  soil  properties  the  process  of  lens  formation  can  alternate  with  frost 
penetration. 

In  sect.  4  we  examine  the  model  when  the  temperatures  are  assigned  at  the  boundaries  of  the  soil. 
Conditions  involving  the  boundary  data  and  the  known  properties  of  the  soil  in  order  to  discriminate 
the  kind  of  process  have  been  obtained  also  in  this  case.  Moreover,  taking  boundary  temperatures 
depending  on  time,  the  phenomenon  of  ’’banding”  (formation  of  several  lenses)  is  correctly  simulated 
by  the  model. 
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1.  The  mathematical  model 

1-1  Basic  assumptions 

Let  us  consider  a  sample  of  soil  saturated  with  water  (in  both  phases)  and  subjected  to 
imposed  gradients  or  temperatures  on  the  top  and  on  the  bottom  of  the  soil. 

We  observe  three  different  regions:  an  unfrozen  part,  from  the  bottom  up  to  the  freezing 

front,  a  frozen  fringe  with  grains  of  soil,  water  and  ice,  and  a  frozen  part,  with  ice  and 
grains. 

In  the  frozen  part  pure  segregated  layers  of  ice  can  form  {ice  lences):  in  this  case,  we 
have  that  the  upper  boundary  of  the  frozen  fringe  (base  of  the  lens)  is  at  rest.  On’ the 
other  hand,  if  a  macroscopic  break  in  the  soil  does  not  occur,  we  can  observe  the 

penetration  of  the  freezing  front  and  the  movement  of  the  frozen  fringe  towards  the  base 
of  the  soil  ( frost  penetration ). 

Many  models  have  been  proposed  in  order  to  describe  the  phenomenon. 

We  are  going  to  give  a  mathematical  description  of  the  freezing  process  following  the 
model  proposed  in  [3]. 

Let  us  take  a  vertical  coordinate  z,  positive  upwards,  and  choose  the  fixed  base  of  the 
soil  as  z  =  0. 

We  consider 

zF{t)  as  the  boundary  between  the  unfrozen  soil  and  the  frozen  fringe; 

*s(<)  as  the  boundary  between  the  frozen  fringe  and  the  frozen  soil; 

zL{t)  as  the  upper  boundary  of  the  forming  (or  just  formed)  ice  layer,  that  is 
immediately  over  the  frozen  fringe 

zt{^)  as  the  top  of  the  soil. 

In  order  to  take  into  account  the  coexistence  of  water  and  ice  in  the  frozen  fringe,  we 

introduce  the  unknown  such  that  eu  is  the  volumetric  ice  content  in  a  unit  volume  of 
mixture  ice-water-soil. 

We  have  i/  =  0in  the  unfrozen  soil,  0  <  u  <  1  in  the  frozen  fringe,  u  =  1  in  the  frozen  soil. 

We  base  our  discussion  on  the  following  statements  (see  [3]  and  [4]  for  more  details): 

(//i)  The  process  can  be  described  by  a  one-dimensional  model. 

(H2)  The  porous  matrix  in  the  unfrozen  soil  is  undeformable  as  well  as  the  skeleton  of 


the  frozen  fringe  (  porosity  e  is  constant  in  0  <  z  <  zs(t)). 

(H3)  In  the  heat  balance  equations  terms  containing  derivative  of  the  temperature  with 
respect  to  time  and  convective  terms  are  neglected. 

(ff4)  Pore  ice  in  the  frozen  fringe  is  at  rest  with  respect  to  the  porous  matrix  (no 
regelation);  moreover,  there  are  not  capillary  interfacial  effects  between  ice  and  water  in 
the  frozen  fringe. 

(ff5)  The  water  flux  qw  is  given  by  the  Darcy’s  law  in  the  unfrozen  soil: 


.  9pw(z,t) 

0  dz 


0  <  z  <  zF{t) 


with  a  constant  hydraulic  conductivity  KQ,  while  in  the  frozen  fringe  the  water  is  driven 
by  a  couple  of  gradients  [5]: 


9 


Si  - 

w 


-KX(T) 


dpw(z,t) 

dz 


-K2(T) 


dT(z,t) 

dz 


zF(t)  <z<  zs(t) 


where  the  empirical  functions  K1  and  K2  are  properties  of  the  soil. 

(H6)  The  water  pressure  pw  at  z  —  zs(t)  attains  in  any  case  the  value  cr  >  0,  that  takes 
into  account  the  overburden  pressure  loaded  over  the  soil.  At  the  base  of  the  soil  the 
water  pressure  is  the  atmospheric  one  and  pw  is  continous  through  the  boundary 

z  =  zsit)' 

( H7 )  The  boundary  zF(t)  is  the  isotherm  T  —  0  {T{zF{t),t)  =  0);  moreover,  the 
temperature  T  is  continuous  with  respect  to  the  spatial  coordinate  z  throughout  the 
whole  soil. 

( H8 )  The  ice  volumetric  content  v  =  i/(T)  (volume  of  ice  in  a  unit  volume  of  mixture)  is 
a  given  decreasing  function  of  the  temperature  T,  such  that  z/(0)  =  0  (this  means  that  at 
z  =  zF(t)  there  is  no  change  of  phase). 

(i/g)  The  quantities  ku  k-  (thermal  conductivities),  pw ,  p p8  (specific  densities)  and 
L  (latent  heat  per  unit  volume)  are  constant.  Suffixes  u,  //,  /,  w ,  «  and  5  refer 
respectively  to  unfrozen  soil,  frozen  fringe,  frozen  soil,  water,  ice  and  soil  grains.  On  the 
other  hand,  the  termal  conductivity  in  the  frozen  fringe  k ^  depends  on  the  temperature 
T.  Empirically,  k ^  is  determined  as  a  function  of  i/,  that  is  in  turn  a  function  of  T  alone 
by  virtue  of  hypothesis  ( H8 ). 

The  set  of  equations  which  describe  the  freezing  process,  valid  both  in  case  of  lens 
formation  and  in  case  of  frost  penetration,  is  the  following: 


-2- 


(1.1) 

(1.2) 

(1.3) 

(1.4) 

(1.5) 


?2)«' 


d{ST(z^ 

tA  dz  1  ~ 0 


qm=-K0^£ 

Tz(kfAT^)  =  o 
«#(*>*)  =  -  Kx{T{z,tyfe*^ 
AfdT^t)\  n 

M  dz 


K2(T(z,t)) 


fgW) 

dz 


with  the  boundary  conditions 


}0  <  z  <  z p(t) 

( unfrozen  soil ) 


>  zF{t)  <z<  zs{t) 
|  ( frozen  fringe ) 

zsW  <  z  <  zx(t) 
(. frozen  soil) 


(1.6) 


0T(O,<) 

-  -«oW  or>  alternatively, 


(1.6a)  T(0,0  =  h(t) 

(L7)  PjO,<)  =  0 

C1-8)  ITM<),<)J  t  =  T(zp(t),t)  =  0 

(1.9)  ! 

(L1°)  -C(<)  =  o 

O-11)  Ip«,(*f(*)><)J±  =o 

(L12)  mzs(t),t)j±=0 

(1.13)  Lpwgw(zs(t)^t)  =  (1  -  Vs)eoPwLzs(t)  +  [  -  k-T(zs(t)ii) j  + 

(1.14)  Pw(zs(t),t)  =  o 


(U5)  PizT(t)  =  Pyj(zs(t),t)  +£0(1  -vS)(Pi-pJzS(t) 
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(1-16)  — ^  ^  =  ~ai(0  or>  alternatively, 

(1.16a)  T{zj<(t)jt)  =  g(t) 

with  the  initial  conditions 

(1.17)  «s(°)  =  fc 

(1.18)  zT(0)-H>b. 

The  symbol  [y] t  in  (1.8),  (1.9),  (1.11)-(1.13)  denotes  the  ’’jump”  of  the  quantity  y 
passing  through  the  height  z: 

Ml  =X+  ~X~  =Hm,x(z)-Hm_x(z) 
z-*z + 

(obviously,  z-*z  +  and  z-*z~  respectively  mean  left-hand  and  right-hand  limit,) 

In  (1.9)  and  (1.13)  the  thermal  conductivity  k  takes  the  appropriate  index  u,  //,  i  or  /. 
We  comment  briefly  the  set  of  equations  (1.1)-(1.18).  We  refer  to  [9]  for  more  details. 
Equations  (1.1),  (1.3),  (1.5)  are  the  heat  balances  in  the  hypotheses  (H3)  and  (// -, 0 ) .  We 
remark  that  the  latent  heat  L  does  not  appear  in  the  heat  balance  (1.3),  since  we 
eliminated  the  derivative  of  the  temperature  T  with  respect  to  time  t.  In  pratice,  we  are 
assuming  that  the  energetic  contribution  due  to  the  change  of  phase  in  the  frozen  fringe 
can  be  neglected,  while  the  most  remarkable  effects  of  latent  heat  release  are  on  the 
boundary  z  =  zs(t),  where  (1.13)  holds. 

Equations  (1.2)  and  (1.4)  come  from  ( Hs ).  Equations  (1.6)  and  (1.16)  prescribe  the 
thermal  fluxes  -a0(<),  —  cq(<)  (where  a0,  cq  >  0)  at  the  base  and  on  the  top  of  the  soil, 
respectively.  Alternatively,  the  boundary  temperatures  h(t)  >  0,  g(t)  <0  at  the 
extremities  of  the  soil  can  be  assigned  (equations  (1.6a)  and  (1.16a)).  The  initial 
conditions  (1.17)  and  (1.18)  give  the  initial  position  of  the  freezing  front  zs  and  the 
initial  height  .of  the  soil,  respectively.  Notice  that  the  initial  thickness  of  the  frozen 
fringe  is  unknown,  since  the  initial  position  of  the  isotherm  zF  is  not  given. 

Equations  (1.7),  (1.11)  and  (1.14)  derive  from  ( H6 ),  while  (1.8),  (1.9)  (heat  balance  at 
zF),  (1.10)  (mass  balance  at  zF )  and  (1.12)  come  from  (HT). 
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Equations  (1.13)  and  (1.15)  axe  respectively  the  heat  and  mass  balances  at  „  (we  put 

The  unknow  quantities  of  the  problem  are  the  temperature  T(z,t),  0  <  z  <  zT(t),  the 
water  pressure  PJz,t),  0  <z<zs(t),  the  hydraulic  fluxes  £(,),  0  <z<zF{t)  and  qfj(z,t), 
zF(i)  <z<  zs(t ),  the  free  boundaries  z  =  zF(t),  z  =  zs(t)  and  z  =  zT(t). 


1.2  Temperature 


Integrating  in  each  region  of  the  soil  (1.1),  (1.3),  (1.5)  and  taking  into  account  of  the 
boundary  conditions  (1.6),  (1.8),  (1.9),  (1.12)  e  (1.16),  or  (1.6a),  (1.8),  (1.9),  (1.12)  and 
(1.16a)  we  get  the  following  formulas  for  the  temperature  T(z,t). 

In  the  case  of  prescription  of  heat  fluxes  at  the  extremities  of  the  soil  (equations  (1.6) 
and  (1.16))  we  have 


-  <*o(f) 


(1.19)  = 


J  *//(  T(z,i)) 
I  kua o(0 


~ai  (0 


From  (1.19)  one  gets: 


0  <  *  <  zF(t) 


zf({)  <z  <  zs{t) 


zs(t)  —  2  ^ 


(1.20) 


-<*omz-zF(t)) 

T(z,t) 

'  J  kff(r})drl  =  ~  kua0 ■(<)(*  -  zF(t)) 
0 


T(z,t)=  —  (<)(z  —  Zs(t))  +  Tg(t) 


^  ^  2  <  zf(*)  ( unfrozen  soil ) 

zF(t)  <z<  zs(t)  (frozen  fringe) 
zs(i)  <  7  <  zT(i)  (frozen  soil) 


where  Ts(t)  =  T(zs(t),t)  satisfies 
Ts(t) 

(L21)  J  kff(r,)dr,=  -kua0(t)(zs(t)-zF(t)). 

0 

On  the  contrary,  if  we  prescribe  the  temperatures  at  the  base  and  on  the  top  of  the  soil 
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(equations  (1.6a)  (1.16a)),  we  have  the  following  equations  for  T(z,t): 


(1.20a)  { 


T(z,t) 


zF{t) 


J  kff(w)dT)  = 


z_£ '(t)  ~  ; 

ZF(t) 


Km 


T(z,t)  ~  *s(t)  -%\t)Ts(t) + -m)9(t) 


where  Ts(t)  is  such  that 


Q<z<  zF(t) 


zF(t)  <z<  zs(t ) 


z s^')  <  z  —  ztW 


( unfrozen  soil ) 
( frozen  fringe ) 
( frozen  soil) 


Ts(t) 

(1.21a)  J  kjj(r))drj  ~  -F---~f-\uh{t). 
o  F 

We  remark  that,  as  a  matter  of  facts,  we  would  have  to  introduce  in  the  frozen  soil 
further  boundaries  separating  the  layers  of  pure  ice  form  the  layers  of  soil  grains 
saturated  with  ice.  In  [9]  we  showed  that,  if  we  assume  that  k{  &  k^  the  ” diffractions”  at 
those  boundaries  can  be  eliminated  and  the  temperature  T  in  the  frozen  soil  is  a  straight 
line  with  respect  to  z,  as  in  (1.20)  and  in  (1.20a). 

1.3  Water  flux  and  pressure 

Integrating  (1.2)  and  using  (1.7)  we  obtain 

(1.22)  Pw(*,t)=-k&*,  0  <z<zF(t). 

On  the  other  hand,  if  we  neglect  in  the  mass  conservation  equation  in  the  frozen  fringe, 
that  is,  under  the  hypotheses  (JJ2)  and  (ff8) 

(1.23)  e0(Pt  -  Pw)§±  +  Pwf/J(z,t)  =  0, 

we  neglect  the  term  with  the  time  derivative,  in  the  spirit  of  the  quasi-steady  approach 
(Uj),  we  get,  taking  also  into  account  (1.10): 

(1.24)  qfJ  =  qfJ(t)  =  q“(t),  zF(t)  <  z  <  zs(t). 
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We  call  qjt)  the  function  that  represents  the  hydraulic  flux  at  each  time  t  from  the 
base  of  the  soil  2  =  0  up  to  the  boundary  z  =  zs(t). 

From  (1.13),  (1.15),  (1.19)  and  (1.24)  we  deduce  the  mass  and  heat  balances  at  the  front 
in  the  following  form,  in  the  case  when  the  thermal  fluxes  e*0  and  a,  are  assigned: 

(L25)  Lpwgw(t)  =  (1  -  vs)epwLzs(t)  +  kfai(t)  -  kua0(t). 

(1.26)  P{zT(t)  =  pwgjt)  +  e0(l  -  vs){pt  -  pw)zs(t) 


In  a  similar  way,  if  the  boundary  temperatures  h  and  g  are  prescribed,  (1.26)  still  holds, 
while  instead  of  (1.25)  we  have: 


( 1-25.)  Lpw,M  =  (1  -  ,s)CfwLism -  tM~Ts W  _  W> 

JzT{i)-zs{t)  zF(t) 

By  means  of  (1.19)  and  (1.4),  one  finds  the  following  formulas  for  the  water  pressure 
gradient  in  the  frozen  fringe: 


(1.27) 


dz 


_  1  ,  v 

Kl  kff(T(z,t))  -9w(*) 


in  the  case  (1.6),  (1.16); 


zf(^)  —  z  <  zs(0 


(1.27a) 


gPu,(M)__  1  ( K2(T(z,t))kuh(t) 

dz  K^T^t))  [  kff(T(z,t))zF(t)  ~ 


Zf(0  5:  Z  <  Zg(t) 


in  the  case  (1.6a),  (1.16a). 

Let  us  assign  now  (1.6)  and  (1.16)  and  integrate  (1.40),  with  the  boundary  condition 
(1.14),  we  get,  making  use  also  of  the  second  equation  of  (1.19): 


(1.28) 


Tf >^)-uoW2 
M  =  j  - ff{r,) 


T(z,t) 


*ifo) 


■dr], 


zF{t)  <z<  zs(t) 


Imposing  the  continuity  of  pw  through  z  =  zF(t)  (condition  (4.1.11))  and  taking 
account  of  (1.8),  we  find: 


into 
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(1.29) 


TJt)  Ts  .  r^K01 0(0 


J _ ^to(^)  _  /^\ 

— - zF\i) 


Analogously,  if  we  prescribe  (1.6a)  and  (4.1.16a),  equations  (1.28)  and  (1.29)  becomes 
respectively: 


(1.28a)  pw(z,t)  =  cr-  J 
T(z,t) 


Ts^K2(t])- 


kuh(t) 


kjf(v) 


KM 


-dr],  zF(t)<z<zs(t) 


(1.29a) 


TM)  is  /_x 

/  *v'1 — Ofir  "(,) 

"  +  J  - ^5) - d”  : 

0 


■tetMO 


Since  the  water  flux  depends  only  on  time,  we  can  achieve  a  formula  for  qw  by 
evaluating  (1.4)  on  z  =  zs(t): 


C1-30)  ?«,(*)  = 


.Ki(TsW)SkMM  +Kl(Ts(t))j0M_ 


or 


(1.30a)  gjt)=  -Kl(Ts(t)fe£gfM  +K2(Ts(t)}---J“h(% 


z  F(t)k  j  ^(T  s{t)) 


1.4  The  mathematical  problem 

We  have  obtained  in  the  case  (1.6)  and  (1.16)  the  set  of  equations  (1.21),  (1.25),  (1.26), 
(1.29)  and  (1.30)  with  the  initial  conditions  (1.17)  and  (1.18);  let  us  write  again,  by  the 
sake  of  convenience,  the  equations: 


-8- 


TS(t) 

zF(t)  -zs(t)  =  j-l—J  kff{v)dr] 

«Jt)  =  (1  -  *3)e*s(i)  +  kJ^Z^ 

LPw 

PizT^)  =  Pw^wi1)  +  e(l-  Vs)(Pi  - Pw)‘zs(i) 

( sfl )  {  Ts(t) 

r  +  J  =  • 

*s(°)  =  b>  zt(°)  =  H>b. 

On  the  other  hand,  when  the  boundary  temperatures  are  assigned  (case  (1.6a)  and 
(1.16a)),  the  set  of  equations  is  given  by  (1.21a),  (1.25a),  (1.26),  (1.29a)  and  (4.1.30a), 
together  with  the  same  initial  conditions: 


(Sfmp) 


TS(t) 

Pi*T(t)  =  Pwqw{t)  +  e(l  -  US)(P{  -  pjzs(t ) 
Ts{t )  K2(ri)-^tlZ,FP\ 


I 


'  +  '  - - *1  -  2f,(«) 


*l(l) 


<m=-ki(Ts(  +Vs(l) 

z5(0)  =  6,  zr(0)  =  H  >  b 


Remark  1.1.  We  will  assume  that  at  the  boundary  zs  holds 
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where  V{  and  Vs  are  the  velocities  of  ice  and  soil  grains,  respectively. 
In  that  case,  the  heave  rate  is  given  by 


kT  =  V,{zs,t)  +  =  Vs(zs,t)  + 


and  the  porosity  e  +  (zs,t )  can  be  calculated  a  posteriori  by  the  formula  (erf.  [9]) 


£(zsS)  + 


f  T  ~  fof 5 

zT-zs  ‘ 


The  five  equations  of  (SyZ)  and  (Stmp)  contain  the  same  six  unknown  quantities  zF(t ), 
2s(0>  Iwi1)-!  Ts(t)  e  Pw( z><)-  Nevertheless,  we  have  not  yet  introduced  a  criterion  in 

order  to  discriminate  the  case  of  lens  formation  from  frost  penetration. 

Before  discussing  such  as  question,  we  need  to  describe  the  propertis  of  the  given 
functions  i /(T),  KX{T),  K2(T)  e  kff(T). 


1.5  The  functions  v{T),  K^T),  K2{T)  e  k^{T)  in  the  frozen  fringe 

As  for  the  volumetric  ice  content,  we  have  to  distinguish  the  two  possibilities: 

i)  v(T)  <1  V  T  <  0 

ii)  3  T<  0  such  that  v{T)  =  1  /orT  <  T. 

We  will  restrict  our  analysis  examinating  only  the  first  case,  remarking  that  the  results 
we  will  get  can  be  easily  formulated  for  the  second  case. 

We  assume  for  the  function  v  to  have  the  same  properties  (cfr.  hypothesis  (H6)  and  [7]): 
(*lo)  0  <  "(T)  <  1  ,  ^Sr~  <  0  ,  KO)  =  1,  Urn  v(T)  =  1 

ul  —  oc 

If  we  extend  the  definition  of  v  to  the  other  regions  of  the  soil  by  setting: 

v{T(z.t))  =  0  for  0  <z<zF(t) 

u(T(z,t))  =  1  for  zs(i)  <z<  zT(t), 
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We  remark  that  *  is  continuous  through  the  boundary  *  =  zF(t),  by  virtue  of  {H6\  while 
it  is  generally  discontinuous  through  *  :=  zs{t). 

The  function  K2{T)  is  positive  and  decreasing  to  zero  for  T  decresing  (see  [8]): 


(^n)  K^T)  >  0  ,  >  0  ,  lim  KX{T)  =  0 


T-+-< 


Furthermore,  we  assume  that  K2(T)  is  such  that  (see  [6]): 
(H12)  K2(T)  >  0  ,  ~fF'>  >  0  ,  lim  _K2(T)  =  0 


Let  us  comment  the  choise  (Hu).  On  the  ground  of  experimental  data,  K2  is  an 
increasing  function  of  T.  A  question  arises  for  the  isotherm  T  =  0,  that  is  the  lower 
boundary  of  the  frozen  fringe:  from  one  hand,  *2  should  achieve  its  maximum  value, 
since  the  temperature  is  decreasing  with  respect  to  z  (see  (1.20)  and  (1.20a));  on  the 
other  hand,  from  t/(0)  =  0  we  deduce  that  only  water  is  present  at  that  boundary 
Actually,  supposing  that  the  effects  of  the  thermal  gradient  on  the  water  flux  are  due  to 
the  simoltaneous  presence  of  water  and  ice,  it  would  be  more  appropriate  to  assume  for 
r2{T)  to  have  a  maximum  for  some  T  <  0,  and  to  be  decreasing  for  T  >  T  up  to 

*2(0)  =  0.  However,  the  temperature  T  is  so  close  to  zero  that  a  crescent  profile  for 
A2(T)  can  be  considered  as  a  good  approximation. 

The  thermal  conductivity  in  the  frozen  fringe  can  be  written  as  (cfr.  [1]): 

kff(T)  -  (1  ~  v{T))ekw  +  +  (1  -  £)ka 

Since  it  is  known  that  (see,  for  istance,  [2])  kt  >  kw,  we  deduce  from  (Hw): 

o-«) 

The  function  kff(T)  has  its  minimum  for  T  =  0,  where  it  assumes  the  value  ku,  which 
represents  the  thermal  conductivity  of  the  unfrozen  soil: 


ku  —  ekw  +  (1  ~~c)ks 


For  T  decreasing  to  tends  to  k,  (thermal  conductivity  of  the  frozen  soi 


1): 
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kf  =  eki-h(l-e)ks. 


1.6  Discriminating  between  lens  formation  and  frost  penetration 


The  basic  assumption  we  made  about  the  water  pressure  is  (1.14):  at  the  boundary  zs 
pw  achieves,  in  any  case,  the  positive  value  a  depending  on  the  kind  of  soil  and  on  the 
overburden  pressure. 

Let  us  assume  that  a  is  a  constant:  this  corresponds,  essentially,  to  neglect  the  variation 
of  weight  of  the  frozen  part  that  leans  on  zs 
Derivating  (1.14)  with  respect  to  time,  one  gets 


(1.32) 


dPw(zs(t)->i) 

dz 


zs(t)  + 


dt 


—  0. 


In  the  spirit  of  the  quasi-steady  approach,  we  assume  that  the  second  term  in  previous 
fromula  is  negligible;  so,  (1.32)  is  approximated  by 

d-33) 

Therefore,  whenever  z5(*)  <  0,  which  is  a  condition  associated  with  the  process  of  frost 
penetration,  it  must  be: 

(1.34)  dPw(zs{i\t)  _  q  frost  penetration . 


On  the  other  hand,  during  the  process  of  lens  formation  the  front  z  —  zs  is  at  rest 

[9]): 


(see 


zs(t)  =  0,  lens  formation. 

Taking  into  account  of  the  initial  condition  (1.17),  we  get 
(1.35)  zs(t)  =  formation. 

We  show  now  that,  during  the  lens  formation,  the  water  pressure  gradient  must  verify 

3't  z  —  z  S ' 
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(1.36) 


dp^*Kn 

~~Fz  ^  °> 


lens  formation. 


indeed,  if  dpw(b,t)/dz<  0  for  some  <  during  the  process  of  lens  formation,  it  must  exist  a 
point  z  in  the  frozen  fringe  so  that  pw(z,t)>a.  This  is  in  contradiction  with  the  fact 

that  we  assign  to  the  condition  Pu,  =  <r  the  property  to  separate  the  soil.  Thus,  (1.36) 
must  hold.  v 

Let  us  now  examinate  the  case  when  both  the  quantities  in  (1.33)  vanish. 

When  z5  =  0  and  zT  =  VJzs.l)  is  not  zero,  we  have  (see  remark  1.1) 

e  +  (zs,t)  =  l 


that  means  that  only  ice  is  present  immediately  over  the  front  z  =  zs. 

On  the  other  hand,  if  it  were  zc(Tl  -  dPw(zs^)it)  .  a.  ,  5  _ 

,  wcie  2S^IJ-  - -zT(t)  =  0  for  some  t ,  we  would  get 

qw(t)  =  0  from  (1.15).  Assume  that  we  prescribe  the  thermal  fluxes  a0  and  ar  From 

(1.4)  we  see  that  this  is  possible  only  if  a0  =  0.  In  that  case,  a,  is  zero,  too  (see  equation 

(1.25))  and  the  temperature  T(z,t)  vanishes  everywhere.  Furthermore,  from  (1.22) 

(1.27)  we  would  have  p„(z,i)  ^  0  for  0  <  z  <  zs(t).  Hence, the  water  pressure  Pw  can  not  bi 

continuous  throughz  =  z5  owing  to(1.14),  unless  <r  =  0.  At  this  point,  it  is  clear  that 

when  both  z5  and  vanish,  zT(i)  can  not  be  zero  for  a  solution  of  (1.1)-(1.18).  A 

similar  argument  holds  in  the  case  of  precscriptions  (1.6a),  (1.16a). 

Let  us  now  make  the  following  important  comment.  From  (1.27)  (or  (4.1.27a))  and  from 
the  properties  of  the  functions  K2(T )  and  t„(T)  (Hn)  and  (4.1.31)  we  deduce  that  the 

water  pressure  gradient  vanish  at  most  in  one  point  in  the  interval  zF(t)  <  z  <  zs(t).  This 
entails  that  whenever 

( 1  371  ^PtAzs(t)ii)  .  . 


holds,  then  the  water  pressure  can  not  exceed  the  value  a  in  any  point.  Thus  (1  37) 
that  is  a  necessary  condition  in  both  the  cases  of  lens  formation  and  frost  penetration 
and  that  requires  the  knowledge  of  quantities  calculated  only  on  the  front  z  =  zs(t ),  is 
also  a  sufficient  condition  in  order  to  guarantee  a  consistent  profile  of  the  pressure  p 
(m  the  sense  that  Pw  <  „)  in  the  whole  frozen  fringe. 


We  remark  that  if  (1.37)  is  satisfies,  not  only  PJz,t)  <  a  in  the  frozen  fringe,  but  p 
non  decreasing  function  of  2-  for  each  fixed  time  t. 


is  a 


On  the  ground  of  the  previous  analysis,  we  introduce  the  following  criterion  C  of 
discriminating  between  lens  formation  and  frost  penetration: 


(C) 


zs(t)  =  0  together  with 

< 

zs(t)  <  0  together  with 


dpw(*SJ) 

dz 


dpw(zsJ) 

dz 


lens  formation 
—  0  o frost  penetration 


Actually,  if  z5(*)  =  0,  then  e  +  (zs,t)  =  1;  moreover,  the  water  pressure  verifies  pw(z,t)  <  a , 
*  e  [*f(*)?*s(*)]*  On  the  other  hand,  if  zs(t)  <  0,  the  water  pressure  gradient  at  z  =  zs  must 
vanish  (see  (1.33)). 


1.7  Further  conditions 


We  have  finally  to  add  the  following  constraints  to  the  solutions  of  (S^)  or  (S*  ),  in 

order  to  eliminate  solutions  which  can  not  be  accepted  from  a  physical  point  of  view. 
Such  as  conditions,  whose  meaning  is  evident,  are,  for  each  time  t: 


1(1.38)  0  <  zF(t)  <  zs{t)  <  zT(t) 

(1.39)  qji)  >  0. 
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2.  Constant  thermal  fluxes  at  the  boundaries 

In  sections  2  and  3  we  will  discuss  the  solvability  of  system  (Sfl).  We  will  investigate  the 
following  points: 


-  thermal  fluxes  given  at  the  boundaries  of  the  soil  as  costant,  in  the  case  of  lens 
formation  (par.  2.1)  and  of  frost  penetration  (par.  2.2); 

-  thermal  fluxes  at  the  boundaries  given  as  functions  of  time  (sect.  3),  with  special 
attention  to  the  transition  processes  from  one  phenomenon  to  the  other  (par.  3.2). 

2.1  Lens  formation;  a0)  constant 


U‘  “S  aSSign  ‘he  ttermal  fluxes  (L6>  “d  (1-16)  by  taking  n0  and  as  constant  and 
consider  the  system  <S„)  in  the  case  of  lens  formation.  Taking  into  occunt  of  (1.35),  the 
set  of  equations  to  be  solved  is: 

Ts(t) 

f(2.i)  zF(i)=£j_y  kff(v)dn+b 


(S//) 


(2.2) 


0 

__  i  -  kua0 

9w  LP~W 


TS(t), 


(2.3)  c+^+/  =» 

(2.4)  =  _  Xl(TsW)^m  +  ma  t  „ 

s  s‘  ttfJsUy,  “  " 

(2.5)  p.i  r(t)  =  p  j 


where 


(2.6)  efa,«,)=  -C*1  =  «e. 


^Pw^ua  0  0 


The  unknown  quantities  in  (S*)  are  the  boundaries  «F(t),  rT«),  the  temperature  r,(t), 

e  water  flux  gJt)  and  the  water  pressure  Pw(zs(t),t):  the  initial  condition  is  (1.18): 
zT(  0)  =  H>b. 
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Moreover,  the  criterion  (C)  introduced  in  par.  1.7  imposes 

(2-7> 

Conditions  (A)  (equations  (1.38)  and  (1.39))  reduce  to 

(2.8)  kfa1-kua0  >  0 

(2.9)  0  <  zF(t)  <  b. 

By  using  (2.1),  we  see  that  (2.9)  is  equivalent  to 

(2.10)  -kua 0b<  J  kff(Tj)dr}<  0. 

0 

In  particular: 

^  =  0  if  and  only  if  /  kff(v)dr)=  -ku<*0b; 


(2.11) 


zF  =  6  if  and  only  if  Ts  =  0 . 


Condition  (2.7)  can  be  written  in  the  following  way: 


(2.12)  ckff(Ts)<K2(Ts),  c>  0. 


In  other  words,  we  have  to  find  Ts  satisfying  the  equation 


(2.13)  G(Ts)  —  0 


where 


(2.14)  6?(s)  =  a  +  c 


V*o] 


5  ^ 

t  +  +/ 


'KM)  -ckff{y) 

KM 


■dr} 


and  such  that  (2.10)  and  (2.12)  are  satisfied. 
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It  is  immediately  seen  that  a  solution  of  (2.13)  is  constant,  since  no  term  depends  on 

,  furthermore,  the  solution  Ts  can  not  be  positive,  otherwise  condition  (2.10)  would 
be  violated. 

Consider  now  the  function  K2(s)-ckff(s),  s  <  0,  c  >  o. 

By  virtue  of  the  properties  of  the  functions  K](.)  and  t„(s)  (see  (ffn)  and  (1.31)  in  par. 

/’  “  ”  e“1,Iy  seen  t,iat  for  ;in'  re(0,A'2(0)/tJ  there  exists  exactly  one  temperature 
Ts(c)  such  that 

(2.15)  K2(T*s(c))  =  ckff(T*s(c)). 

Moreover,  we  have: 


(2.16) 


Hm  T*s(c)  =  -  oo 
c-f  0 


lim  t*s{c)  =  0 ■ 


Owing  to  (2.12),  the  solution  Ts  of  (2.13)  has  to  be  searched  in  the  interval  \T'M)  o] 

wth  c  e(0,n2(0)/lu).  The  lenght  of  that  interval  depends  obviously  on  the  data  «  !' 

through  c.  Furthermore:  07  1 

— ^ —  =  0  if  and  only  if  Ts  =  Tj(c). 

If  C  >  K2(0)/ku,  there  are  not  solutions  of  (2.13)  so  that  (2.12)  is  fulfilled. 

On  the  Other  hand,  condition  (2.10)  means  that  the  solution  Ts  of  (2.13)  must  be  in  the 
interval  [r5(<*0),oj,  where  Tbs(a0)  is  the  one  value  such  that 
Tbs 

(2.17)  -f~ua0b=  J 
0 

Since  fkf/(T])dT)  is  an  increasing  function  of  s,  it  is  easy  to  verify  that  to  every  <*0  there 
corresponds  a  unique  Tbs(a0 );  moreover: 

(2J8)  a^m+00T^o)=—  • 

We  remark  that  zF  =  0  if  and  only  if  Ts  =  Tbs(a0). 
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In  order  to  discuss  the  solvability  of  (2.13),  we  must  consider  the  two  possibilities: 


0  T*s(c K,«i))  >  Ts(«0) 

*"*')  rs(c(«o,ai))  <  Tg(a0) 

We  call  (I)  the  set  of  equations  (2.10),  (2.12)  and  (2.13)  and  we  state  the  following 
LEMMA  2.1.  System  (T)  has  exactly  one  solution  Ts  if  and  only  if 

(2.19)  Cr(T^(c(a0,a1))  <  0  in  the  case  i), 

(2.20)  G(Tbs(a0))  <  0  in  the  case  ii). 


Moreover ,  Ts  is  constant . 

Dim. 

By  derivating  (2.14)  with  respect  to  s  one  finds: 


G'(s)  = 


(Ko(K2  ~  ckff(s))  +  ckff(s)Kl(s) 

Kfs) 


We  see  that  G'(s)  >  0  if:r£(c)  <  s  <  0.  Furthermore,  G(0)  =  a  +  c  0.  The  solution 

Ts  of  (1)  must  lie  in  the  interval  [T^,0],  where  =  T^a^af)  =  ma  dTs(“o)»nw«o, «!))}■ 
We  deduce  that  there  exists  a  unique  solution  of  (T)  if  and  only  if  G(T <  0,  that  is 
equivalent  to  (2.19)  and  (2.20). 

Notice  that  in  (2.19)  equality  holds  (in  that  case  the  solution  of  (2.13)  is 
Ts  =  T5(c(a0,a1)))  if  and  only  if  the  water  pressure  gradient  computed  at  z  =  b  vanishes. 

In  (2.20)  equality  holds  if  and  only  if  zF  =  0  (the  base  of  the  soil  and  the  lower  boundary 
of  the  frozen  fringe  are  the  same  line).  In  that  case,  the  solution  of  (T)  is  obviously 
Ts  =  Tbs(a0).  □ 


Our  next  goal  is  to  locate  the  points  (a^aj)  on  the  quarter  of  plane  {(a0,ai):  a0  >  O.cq  >  0} 
such  that  condition  (2.19)  or  (2.20)  of  the  previous  lemma  is  verified.  The  region  that 
we  will  single  out  and  that  we  call  £,  will  correspond  to  all  the  pairs  of  thermal  fluxes 
(oojQfj)  and  only  those  such  that  the  process  of  lens  formation  is  induced. 

Since  ce(0,K2(0)/ku),  the  region  i.  is  contained  in  the  angle  /?: 
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(2.21)  l  c/3-  j(<*0,a1):  a0  >  0,  kua0  <  kfar  <  kuaj^l  +^k2( 0))) 

We  start  by  remarking  that  the  distinction  between  case  i)  and  case  ii)  is  nothing  but  a 
partition  m  the  angle  p.  Indeed,  let  us  consider  the  set  T*(c(a0,ai))  =  Tbs(a0),  which  we 
can  write,  owing  to  (2.17),  as 

Ts(c(a0>al)) 

(2.22)  (C)  -kua0b=  J  kjf(r])dr) 

0 

easy  to  check,  by  using  (2.16),  that  the  set  (C)  is  a  curve  dividing  p  into  two  parts, 
corresponding  to  the  sets  T*s(c(a0,ai ))  £  Tbs(a0). 

As  a  matter  af  fact,  let  us  take  the  half-straight  lines  starting  from  the  origin 

(2.23)  rc  =  |(a0,ax):  a0  >  0,  kfa1  =  kuaQ(  1  +  pwLc)j, 

with  0  <  c  <  K2(0)/ku.  Going  along  one  of  those  lines  and  letting  a0  grow  from  .zero  to 
+  oo,  we  see  that  the  quantity  T*s  is  a  constant  negative  value,  while  Tbs(a0)  decreases 
monotonically  from  down  to  -oo  (cfr.  (2.17)).  Hence,  there  exists  a  unique  value  a0 
such  that  the  equality  T*s(c)  =  Tbs(a0)  holds;  in  other  words,  there  is  a  unique  point  of 
intersection  P(c)  between  the  straight  line  corresponding  to  c  and  the  set  C.  By  (2.16) 
and  (2.18)  we  have 

(2.24)  |P(c)-0|  =  oo,  Km  |P(c)-O|=0 

"(tt) 

where  the  point  O  is  the  origin  of  axes  in  the  (a0,a1)-plane. 

Moreover,  since  T*s(c)  is  an  increasing  function  of  c,  we  see  that,  as  c  increases  from  zero 

up  to  K2( 0)/ku,  the  distance \P(c)-0\ decreases  monotonically  from  oo  to  zero. 

Let  us  denote  now  by 

(2.25)  =  {(Ofcflfj)  e  0: :  T‘(c(a0,ai))  <  Tbs(a0)J 

(2.26)  %2  ={(a0,ai)  e  p:  T*s(c(a0,ai))  >  T*  C^)} 
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the  two  parts  separated  by  the  curve  C.  The  region  9^  is  bounded  by  the  straight  line 
kfai~kuao  fr°m  below  and  by  the  curve  C  from  above,  while  %2  is  bounded  by  the 
curve  C  from  below  and  by  the  straight  line  k^al  =  kua0(l  +  pwLK2(0)/ku )  from  above. 

By  virtue  of  lemma  2.1  we  can  say  that  the  process  of  lens  formation  occurs  if  and  only 
if 


Tbs 


(2.27) 


<T  + 


/ 


K2{t))  -ckff(t]) 

Ki(l) 


dr\  <  0 


whenever  (a0,<*j)  €  9^, 

TUc) 

(2.28)  .  +  ™>(rJW+^)+/ 


whenever  (a0,ax)  e 
If  we  call 


i  ~  {(ao?ai)  e  (4.2.27)  is  verified} 

^  (4.2.28)  is  verified} 

we  may  conclude  that  JL  =  ^  ^  is  the  part  of  the  angle  j3  where  lens  formation 

occurs. 

We  begin  with  the  following  result. 

LEMMA  2.2.  Let  Ta<  0  be  such  that 
<2-29) 

0 

then  there  exists  a  unique  c,e(0;  K2(0)/kJ  with  the  property  that 
T*s(c, ) 

(2-30)  a+j  - - </,  =  0. 

0 

Dim.  We  first  remark  Jhat  the  function 

f^=cr+J  t$, \d” 

0 
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is  strictly  increasing  and  /,(0)  >  0.  Thus,  the  assumption  (2.29)  is  equivalent  to 
(2.31)  Jjrn  oo/1(s)  =  £  <  0. 

Indeed,  if  (2.31)  holds,  there  is  a  unique  Ta  such  that  (2.29)  holds  and  viceversa. 

Let  us  now  consider  the  family  of  functions  depending  on  the  parameter  c: 

/c(s) "  ^  +  Vs  ^  °>  0  <  c  <  K2(0 )jku. 

0 

We  have: 

(2.32)  / c(s)  >  /j( 5),  Jim+  fc(s)  =  /j(s)  Vs  <  0  (punctually). 

Moreover,  according  to  the  properties  of  h\  and  kff  and  since 

dfM  K2(s)  -ckJf{s) 

&s  Kx(s)  ' 5 

we  see  that  for  any  e  fixed  in  [0,tfs(0)/ij,  the  function  achieves  its  minimum  value 
for  «  =  T*s(c),  where  T*s(c)  is  defined  by  (2.15): 

(2-33)  fc(s)>fc(T*s(c))  Vs<0,  ce[0,K2(0)/ku]. 

Suppose,  contrary  to  our  claims,  that  /C(TJ(C))  >  0  Vc6(0,  K2(0)/ku]  and  consider  the 
case  £  >  —  oo. 

By  virtue  of  (2.31),  once  we  have  fixed  0  <  e  <  - if  2  we  could  find  s  <  0  such  that 

h  + 

On  the  other  hand,  by  (2.32)  we  could  find  c  small  enough  that 
0  </s(S)-/1(S)<e. 

Therefore: 

/?(«)  <  £  +  2e  <  0. 
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On  the  contrary,  from  (2.33)  we  deduce: 


4(*)  >  >  o, 

and  we  obtail  a  contradiction.  The  case  ^=-00  can  be  proved  with  slight  modifications. 
Thus,  we  may  say  that  there  exists  c  e  (0,  K2(§)/ku)  such  that  /-  <  0. 

Since  the  function 

T*(c) 

r  K2(rl)  , 

B(c)  =  a+ j  - - dlh  c  e  (0,  A'2(0)AJ 

0  1 

is  continuous  and  verifies 

T*s(c) 

B(K2( 0)/*J  =  00,  B(c)  <  0,  B'(c)  =  -  J  dr 1  >  0  per  c  e  (0,  K2( 0)/ku], 

o  1 

we  conlcude  that  there  exists  a  unique  cjcce(0,  K2(0 )/ku)  with  the  property  B(c{)  —  0, 
and  the  lemma  is  proved.  □ 

The  following  result  locates  the  region  JL  on  the  (a0,aj)-plane. 

PROPOSITION  2.1. 

i)If 

(2.34)  c  +  J  >  0  Vt?  e  ( —  oo,0)  , 

0  1 

then  JL  =  0. 


ii)  If  3  T q.  <  0  such  that 

T<r 

(2-29)  ’+!  Wfn=^ 

0 

then  JL  is  the  non-empty  set  bounded  by 

(2.35)  dJL  =  |(a0,a1)  :  —  kua0  ,  a0  >  a0  j  U  C1  U  C2 
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where  S0 

is  defined  by 

Ta 

(2.36) 

_  b  h  _  f  ^2^)  i 

w-y  km  * 
0 

and 

TbS(a  0) 

(2.37) 

Cl={(a0>al):  <7+  / 

J 

0 

(2.38) 

C2  ~{(a0)al):  a  +  c  xj 

~  0,0  <  c  < 

K2{q)  ~cKJr})  ^ 

Y^)  dt)  =  0,0  <  c  <  C/j 


Ts(c )  N  Ts(c) 

*H+/ 


-dr]  — 


are  two  curves  contained  in  and  %2  respectively;  in  (2.37)  and  in  (2.38)  c,  is  the 
value  defined  in  lemma  2.1  and  verifying  (2.30). 


Dim.  Case  i)  is  immediately  proved:  if  (2.34)  holds,  equation  (2.13)  can  not  have  non¬ 
positive  solutions  and  the  condition  (2.10)  is  violated.  Hence  2  =  0. 

Let  us  now  pass  to  case  ii). 

We  consider  once  again  the  straight  lines  (2.23).  We  know  that  there  exists  a  unique 
point  P(c)  whose  coordinates  (a0,aj)  axe 

Ts(c(s< 0>«l)) 

(2.39)  -kuaQb~  J  kff(r))di 7,  =  £u50(l  +  pwLc)  ,  per  0  <  c  <  — 

0  k“ 


Conventionally,  we  can  assume  that  5o(0)  =  5,(0)  =  +oo,  recalling  that  the  straight  line 
ro  :  kf<x  1  =  kuao  is  an  asymptote  for  the  curve  C. 

We  exanimate  first  the  region  Evaluating  the  quantity  on  the  left-hand  side  in  the 
inequality  (2.27)  along  one  of  the  lines  rc,  where  c  is  constant,  we  get: 

TsK), 


(2.40)  f(« „,*)  =  „+/  MLzfnM 

J  Ki(v) 


-drj 


(«0>al)  6  ,  0  <  C  < 


By  virtue  of  assumption  (2.29),  lemma  2.2  and  remarking  that  £ F(a0(c),c )  >  0,  we  have 
(2'41)  lJ^oFM=  <^>0  Vce  [0,tf2(0)/*u] 
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(2.42)  F(a0(c),c) 


1: 


0  if  0  <  c  <  C[ 
0  if  c  —  ct 


K2(°) 

>0  ifc,<c<  ' -■ 


Taking  into  account  of  (2.41),  (2.42)  and  noticing  that  F(a0, 0)  decreases  along  the  lines 
rc,  that  is 


(2.43)  dZM<0,  (ao.ariJe*! 

we  conclude  that  whenever  (a0,aq)  e  3^,  if  c{  <  c  <  K2(0)/ku,  there  are  not  any  points 
(o^cq)  on  the  straight  line  rc  such  that  (2.27)  is  verified.  In  other  words: 

(2.44)  %i^f1rc  =0,  Ci<c<  — . 

On  the  other  hand,  when  0  <  c  <  cb  we  have  that  (2.27)  evaluated  on  each  line  rc  is 
valid  if  a0  is  such  that 


(2.45)  c*0  c  K(c),or0(c))  , 

where  ^0(c)  is  the  value  (which  exists  unique  owing  to  the  continuity  of  the  function  F 
with  respect  to  a0  and  to  the  formulas  (2.42))  with  the  property 


(2.46)  ^(aofcjjc)  —  0,  0  <  c  <  cj 

This  is  equivalent  to  state  that  for  0  <  c  <  cl  we  have 

(2.47)  e  n  rc  =  {(ao-o-j):  50(c)  <  a0  <  a0(c)  ,  kfat  =  kua0(  1  +  pwLc )} 
When  c  =  0,  it  is  ao(0)  =  H-oo,  50(c)  =  a0  ,  where  S0  satisfies  (2.36).  Thus: 

\«nr0  =  {(Vl):  a0>SO:  kfa l  =  ^uao}- 


Obviously 
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„<U^N-,e. 

The  set  of  points  in  3^  satisfying  (2.46),  that  is 

{(a0>Qi)  e  a0  -  a 0(c)  ,  kfax  =  fcu50(c)Aua0(l  +  pwLc),  0  <  c  <  c,|, 

make  a  curve  Cx  that  is  the  boundary  of  the  set  «, 

The  curve  Cx  matches  C  defined  by  (2.22)  only  at  the  point  P  =  (&0,&x)  with  coordinates 

T*s(c,) 

(2.48)  ~kuaQb=  J  kff(r,)dr],  kfa1  =  kua0((l  +  pwLc,) 

0 

and  the  straight  line  kf ox  =  kua0  at  a0  -  S0,  where  30  is  defined  by  (2.36).  Moreover,  the 
tangent  to  the  curve  Cx  at  P  alia  curva  Cx  is  the  straight  line 

rc,  ~  {(Q0)al):  ao  >  */ai  =  ku<x0(l  +  pwLc{)j. 

We  will  examinate  now  the  region  %2.  In  order  to  find  the  points  (a0,<*x)  in  %2  which 

have  the  property  (2.28),  we  check  the  values  achieved  by  the  auxiliary  function  H ,  that 

we  are  going  to  define,  along  each  of  the  straight  lines  rc,  0  <  c  <  K2(0)/ku,  where  Tl(c)  is 
constant: 

:V  f  ^2(??)-^//»?). 

/+  J  #,(77)  dri'  (Q0ial)6%2 

/  n 


(2.49)  H(a0,c)  =  *  +  MEs^y>( TUc)  +  r  'KM-Ckff{ri) 


TTI'J 


0 

We  can  exclude  the  value  c  =  0,  since  the  straight  line  r0  :  kfa,  =  Vo  verifies  %2  n  r0  =  0. 
We  have  to  find  the  part  of  %2  where  G(a0.c)  <  0,  owing  to  (2.28). 

It  is  immediately  seen  that: 

(2.50)  H(a0 ,c)  =  F(a0(c),c)  0  <  c  <  K(0)/ku. 

moreover,  along  one  of  the  line  r  we  have: 


(2.51) 


lim  G(aa.c)  =  +00  0  <c  <  K  <0)/k 

(Q0>«l)erc,a0-»  +  oo  ^  C  -  A2W/V 

Arguing  as  in  the  case  i),  we  conclude  that  whenever  (a0,ax)  €  %2: 
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(2.52) 


if  ct<c  <  K2/kuJ  then  ff(aQ)c)  >  0; 
<  if  c  =  clf  then  H(ocQ)c)  —  0; 

if  0  <  c  <  c;  then  H(a0)c)  <  0  when 


(2.53)  a0  €  [  <*0(c),ar0(c)], 


where  a0(c)  has  the  property 


(2.54)  tf(a0(c),c)  =  0,  0  <  c  <  cz. 

The  existence  and  the  uniqueness  of  s0(c),  0  <  c  <  c{  come  form  the  continuity  of  the 
function  H  and  from  (2.51)  and  (2.52).  Thus,  in  %2  we  have: 

\  l  n  rc  =  0,  if  Cl<c<K2/ku 

%2,£nrc  ={(a05ai):  50(c)  <  a0  <  5 0(c)  ,  k jOc^  =  kua0(  1  +  />„,!£:)  j  ,  if  0  <  c  <  cr 
Equation  (2.54)  defines  a  curve  (C2) 

{(a0:al)  e  %2:  <*0  =  5o(c)  )  kja  1  =  *  J*0(c)(l  +  PwLc)->  0  <  C  <  C,| 
that  is  the  one  we  defined  in  (2.38).  □ 


Remark  2.1.  The  point  P  (see  (2.48))  lies  on  the  curve  C2.  Furthermore: 


(2.55)  a0(c<)  =  a0(cj)  =  a0(c,)  =  a0. 

(2.56)  lim 

<*o(c)  =  +  °°- 

c-4)  + 

If  we  denote  by  P,  Px  and  P2  the  points 


P{c)  =  |a0(c),  ^  a0(c)(l  +  PwLc)} 
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( 


(2.57) 


for  each 
increasse 
we  have: 


<  Px (c)  =  |a0(c),  ^  50(c)(l  +  pwLc)j 
P2(c)  =  |50(c),  I0(C)(l+/vtc)J 

c  in  the  range  0<c<C/  (conventionallly,  P2( 0)  =  +oo),  we  see  that  \Px{c)-0\ 
if  c  increases,  while  \P(c)-0[  \P2{c)-0\  decrease  if  c  increases,  besides  that, 


(2.58)  |  Px{c)  -  O  |  <  j  P(c)  -  O  |  <  |  P2(c )  -  O  | 


and  equality  holds  if  and  only  if  c  = 

The  region  L  is  drawn  in  fig.  4.1. 
Remark  2.2.  When  the  point  (a^cq)  ec1? 


the  solutions  of  (S^)  have  the  properties 


(2-59)  zF  =  o,  rs  =  r|(a0), 


(“o^ileCj. 


On  the  contrary,  when  (<*0,a1)6e2  we  find  solutions  of  (sj{)  such  that 


(2.60) 


TS  ~  T5(c(a0,aa)), 


dpw(b,t) 

dz 


(a05tti)  c  e2. 


Once  Ts  has  been  calculated  by  means  of  system  (T),  whenever  (*<,,«,)  ei.  #  0,  We  can 
compute  the  temperature  in  any  point  of  the  soil  by  using  the  formulas  (1.20).  We 
notice  that  the  temperature  T  depends  only  on  the  spatial  coordinate  z  but  not  on  time 

By  (2.1)  we  find  the  (constant)  thickness  of  the  frozen  fringe: 

Ts 

b~Zp  =  k^J  kfMdll 
0 

while  (2.2)  and  (2.5)  allow  us  to  calculate  the  upper  boundary  of  the  soil: 


p  kfa^~  ktlan 


Formulas  (1.22)  (1.28)  give  the  water  pressure  in  the  unfrozen  soil  and  in  the  frozen 
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fringe,  respectively: 


0<*<  zF(t) 


zf( 0  <  *  <  *s(0- 


The  water  pressure  pw  achieves  the  minimum  value  for  z~zF.  Moreover,  owing  to  (2.7) 
and  to  the  properties  of  K2  and  kjp  the  pressure  pw  increases  with  respect  to  z  in  the 
frozen  fringe. 

Recalling  (1.2)  and  (1.27),  we  see  that  the  z-derivative  of  pw  has  a  discontinuity  at 

z  —  Zp\ 


=  i-W  °K- 
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We  investigate  now  the  system  (S/;)  imposing  the  conditions  that  are  peculiar  of  the 
process  of  frost  penetration.  We  keep  the  assumption  a0,  a1  constant.  Recalling  the 
criterion(C)  (see  par.  1.7),  we  see  that  the  set  of  equations  (Sy;)  defined  in  par.  1.5  takes 
the  form: 


TS(t) 


I"  (2.61)  zF(t)-zs(t)=-j±^J  kff(r])dT) 

(2-62)  qjt)  =  (1  -  us)ezs(t)  +  --Lp - 

Ts(t) 

(2.64)  qjt] )  =  JJpj'fyj  kua0 

(2.65)  =  Pwqjt)  +  e(l  -  l/s)(Pi  ~  Pw^sW 


with  the  initial  conditions  (1.17),  (1.18): 


z5(0)  =  b,  zr(0)  =  ff>b. 

The  solution  (Ts(t),zF(t),zs(t),qw(t))  of  (Sfy  must  fulfil  the  following  conditions  (cfr. 
criterion  (C)  and  conditions  (ft)  in  parr.  1.7,  1.8): 


(2.66) 

zs<0 

(2.67) 

0  <  zF(t)  <  zs(t)  <  zF{t) 

(2.68) 

qjt)  >  o. 

Equation  (2.64)  gives  the  water  flux  in  terms  of  the  temperature  Ts(t)  alone.  Expressing 
also  zF  also  zs  in  terms  of  Ts,  one  gets: 


(2.69) 


_  A'0l'yy(T5(<)) 

Zp{)~ 
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(2.70)  zs(t)=- 


^ff^swu  Tfs(t) 

kua0^2(Ts(i))\  J 

'  n 


TsW 

~K^0J  k^)dT]- 

0 

From  (2.62)  and  (2.64)  we  find 


KM 


■dr]  I- 


(2.71) 


h(t)  =  -rkt‘a° 


(i-^ 


MMR  \ 

kfATsW)  ~c) 


where  c(a0,ai)  = 


*/ai  ~  K*o 
LPwkua  o 


Condition  (2.66)  imposes  that  the  solution  Ts(t)  must  verify  at  any  time  t: 
(2.72)  K2{Ts{t))  <  ckff(Ts(t)). 


Differentiating  with  respect  to  time  (2.70): 


(2.73) 


^5(0  = 


_«o./  d  (kfATs(t))\( 
kuao  A>T^K2(Ts(t))J[J  + 


Ts(i ) 

/  kfATs(t)) 

1  kM0’>)+— Kp- 


and  combining  (2.71)  with  (2.73) 
for  the  temperature  Ts(t): 


we  obtain  the  following  ordinary  differential  equation 


(2.74) 

where 


— ~ — ^4^-  <p(Ts(t))Ts(t )  - 
(Kao)  S{  ’  kff(Ts(t)) 


(2'75) 


cr-f 


5 

[K 2 

y  ^ 


*2fo),  ),*//(«) 


/  ; 

We  assume  that  the  given  functions  *„  and  e  ace  regular  enough  (say 

(  ~.°J)  <>rfe  ‘o  guarantee  local  existence  and  uniqueness  of  a  solution  TM  of 

(2.74),  supposing,  for  the  moment,  to  know  the  initial  value  Ts(0).  The  integration  with 
respect  to  time  of  (2.74)  yelds: 
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(2.76) 


f  Ko(1^i's(y))kfj(y)€ 

t  (o)  ^“a°^2^2 (y)~ckff(y)) 


<p{y)dy  =  t. 


Formally,  we  have  found  the  solution  of  (S %),  since,  once  (2.76)  has  been  inverted,  the 
boundaries  zF(t)  and  zs(t)  are  achieved  by  (2.69)  and  (2.70),  respectively;  the  flux  qw(t) 
and  the  surface  zF(t)  are  calculated  by  means  of  (2.64)  and  (2.65),  respectively. 

However,  we  have  to  check  in  which  cases,  that  is,  for  which  values  of  (a0,a,),  the 
solution  is  consistent  with  the  conditions  (2.66)-(2.68).  For  this  purpose,  we  set 

(2.77)  7  =  {(o0)tti):  3!  solution  of  { sf,)  such  that  (2.66)-(2.68)  hold} 

and  we  state  the  following 

PROPOSITION  2.2. 


i)  If  (2.34)  holds,  then  7  =  0. 

it)  If  (2.29)  is  valid,  then  7  is  the  set  bounded  by 

(2.78)  57=  {(a0,a1):a0  =  S0,  aj 

where  SQ  and  ax  has  been  defined  by  (2.48)  and  the  curv  C2  by  (2.38). 
In  particular ,  if  we  define  the  angles 

(2.79)  0O  =  {(a0,aj):  k <  kua0  <  fcua0(l  +  Lpwc{ )} 

(2.80)  /31  =  {(aQ,^):  k^a1  >  £ua0(l  +  Lpwc{)}, 
we  have: 


u)a  if  (a0,a1)e70  =  7n/?0,  then  there  exists  a  time  tf  (finite  or  infinite )  such  that 

hit)  =  0/ 

ii)b  if  (aQ,^)  e  7j  =  7  n  f)v  then  there  exists  a  finite  time  t  such  that  zF(t )  =  0. 
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Dim . 

We  begin  by  evaluating  (2.69)  and  (2.70)  when  t  =  0: 


(2.81)  zF( 0)=  _KoLff(Ts(°))  (  r 

kuaoK2(Ts(0))  +  j 
x  0 


Ts<0) 


km 


dr) 


(2.82)  *s(0)  =  6  =  -  °  ^(r-?(°))  L,  [ 

kuaoK2(Ts(0))[a+ J 
x  0 


TM 


KM 


dv  - 


Ts(0) 


■)  If  (2  34)  holds,  no  temperature  Ts(0)  can  be  founded  so  that  (2.69)  could  be  satisfied. 

'  “  Ts(0>  has  t0  be  non-negative  owing  to  (2.61)  and  (2.67);  but  by  (2.81)  and 

(2.34)  ,t  would  result  (0)  <  0,  in  contradiction  with  (2.67).  Case  i)  is  thus  proved. 

ii)  In  order  to  find  solutions  of  (Sf,)  consistent  with  the  constraints,  it  is  necessary  to 
assume  that  the  temperature  T,  verifying  (2.29)  exists;  by  virtue  of  lemma  2.2 
ere  ex,sts  c,  e  (0,2fj(0)/i J  satisfying  (2.30)  (we  can  consider  the  open  interval,  since  if 

dXt’  °e“SS“ily  C,<KmlK)'  S°'  WC  ^  ‘ha‘  ‘he  “*><*  P-79)  P-80)  »e  well 

According  to  the  properties  of  K,  and  (see  assumption  (if,,)  and  (1.31)),  one  has 


T,  <  T  . 


(2.83)  -jgo  .9  (kff(Ts)\(  7SK2{ij)  \ 

Ts  kua0  dT^K2(Ts)lr+  j  KfifV<  0,  for 

x  0  / 

Since  rF(TJ(e,))  =  o  (cfr.  (2.15)  and  (2.30)),  from  (2.83)  we  deduce  that  the  condition 
ZF  >  0  is  satisfied  if  and  only  if  Ts  <  t£(cj). 

Let  us  now  examinate  the  two  regions  (2.79)  and  (2.80)  separately. 
n)a  (ao,ai)e0o. 


We  remark  that  whenever  (o'0,Qt1 )  e  /?0  the  value  defined  immediately  after 


the  equation 
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(2.71)  fixes  uni  vocally  the  temperature  T*s(c )  such  that  (2.15)  is  verified. 

We  are  going  now  to  discuss  the  solvability  of  (2.82),  that  allows  us  to  achieve  the 
initial  value  of  the  freezing  temperature  Ts(0),  on  the  ground  of  the  knowledge  of  the 
really  assigned  initial  condition  zs( 0)  —  6. 

From  (2.71)  we  deduce  that  zs(t)  <  0  if  and  only  if  Ts(t)  <  Tg(c),  since  f(s )  =  K2(s)/kJf(s ) 
is  an  increasing  function  of  s  and  f(T*s(c))  =  c  (c  is  fixed  once  a0  and  a1  are  prescribed). 
Therefore,  the  initial  temperature  must  verify 

(2.84)  t5(0)  <  T*s(c). 


On  the  other  hand,  because  of  (2.61)  and  (2.67)  evaluated  at  t  =  0  it  must  be: 

T5(0) 

(2.85)  kffWn. 

0 

Condition  (2.85)  is  equivalent  to 

(2.86)  rs(o)  >  Tbs(a0) 


by  virtue  of  (2.17). 

From  (2.84)  and  (2.86)  we  achieve  the  condition 
(2.87)  Tbs(a0 )  <  Tg(c). 


Notice  that  (2.87)  excludes  the  region  3^  defined  by  (2.25): 


,$n%1  =  0. 


Call 

(2.88) 


B(s)  =  - 


kua 0K2(s) 


(  f  KJ',]  i  ,  ’ 

{  J  KM  kiMd^ 


Equation  (2.82)  is 
(2.89)  B(TS(  0))  =  b. 
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It  is  easily  seen  that 


(2.90)  B(T}(c))  >  0,  5(T|(o0))  >  b,  B'(s)  <  0  for  s  <  T£(c). 

We  may  conclude  that  (2.82)  has  a  unique  solution  T5(0)  in  the  interval  fr6„f&n)  T%(c) 
provided  that  S  ' 

(2.91)  B(T*s(c))  <  b. 

It  is  evident  that  B(T*s(c))  =  b  defines  the  curve  C2  (see  (2.38)).  The  region  outlined  by 

(2.91)  corresponds  to  the  right-hand  part  with  respect  to  the  curve  C2  on  the  (a  a  )- 

plane.  05  1 

Remark  2.3.  condition  (2.87)  is  surely  fulfilled  if  (2.91)  holds,  since  the  curve  C  (defined 
by  Ts(a0)  =  Tg(c))  is  closer  to  the  origin  than  C2: 

(2.92)  (K>ai)  €  0o:B(T*(c))  <b}c  {(«„,<*,)  e  0O:  Tbs(a0)  <  T}(c)}. 

Examinating  conditions  (2.66)-(2.68),  we  see  that  the  imposed  conditions  are  fulfilled  for 

<  =  °’  Provided  ttat  (2.91)  holds.  Moreover,  zF{  0)  =  o  if  and  only  if 
r5(0)  =  T‘(ob)  =  Tj(cl). 

Once  the  initial  condition  Ts{ 0)  has  been  calculated  by  means  of  (2.89),  the  temperature 
Ts(t)  is  achieved  by  integration  from  (2.76). 

The  function  <p  defined  by  (2.75)  is  positive  for  sKT,  (cfr.  (4.2.83)).  Furthermore,  since 
r5(c)  <  T*s(ci)<  T„  whenever  (a0,«i)^0,  from  (2.76)  and  from  (2.84)  we  deduce  that 
t5(0)  >  0,  so  fs(i)  keeps  positive  provided  that  Ts(t)  <  T*s(c).  Moreover: 


l^5(t)|< 

It  is  easy  to  check  that  Ts(t)  can  not  have  an  asymptote  in  the  interval  (Ts(0):T*s(c)). 
Thus,  the  temperature  Ts(t)  necesssarily  reaches  by  increasing  the  value  T*s(c)  in  a  finite 
or  infinite  time  tf.  In  the  interval  [o,^)  conditions  (2.66)-(2.68)  are  verified  (see  (2  69) 
(2.70),  (2.72),  (2.73)  and  (2.83)).  When  t  =  tf,  we  have: 


c  W)) 

kff(rsm ) ( u°o) 

kJl-v(Ts(Q))) 
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zp(tf)  -  zHaoiai)  =  - 


Ko  kff(Ts(c)) 


kua0K2(T*s(c)f 


T*s(c) 

V+I  ~ 


\  o 


K\  (v) 


dr} 


i 


zs(tf)  -  zs(a0^l)  -  - 


K0kjf(T*s{c)) 


ku<*0  K2(T*s(c)y 


Ts(c)  K2(t))  ■ 


c  + 


\  0 


I 


*2  (T*s(c)) 

kJf(^s(c)) 


kff(v) 


KM 


dr} 


+ 


T*s(c) 

+^k/  kffWdi>zF(tf) 


vj  })  Ka0c  -  k2(T*s(c))  kuQ° >  0 

where  2s(c(«ojai))  stationary  solution  of  (2.74)  (that  can  not  be  accepted  owing  to 

the  fact  that  condition  (2.66)  is  not  valid),  with  Ts(t)  =  T*s(c)-  We  see  that  zs(tf)  =  0, 
that  is  when  t  =  tf  condition  (2.63)  fails.  By  (2.30)  we  have  zF(t})  >  0  and  equality  holds 
if  and  only  if  c  =  ct. 

Remark  2.4.  The  function  in  the  integral  in  (2.76)  is  positive  for  Ts(0)  <  Ts(t)  <  TJ(c) 
and  the  denominator  vanishes  if  and  only  if  y  =  T*s(c).  Thus,  the  two  possibilities  if  <  oo, 
if  =  oo  occur  according  merely  to 


Ts(c) 

(2.93) 

I 

TS(  0) 

T*s(c) 

(2.94) 

I 

Ts(0) 

K0(l-us(y))kff(y)e  _ 

0 K«o? K2{y)-ckff(y )  ^  y-°° 


kq (1~i/s(y))kjf(y)£  , , , 

— - 2  - - — (fi(y)dy  <  oo 

C kuao )  k2 (y)~ckff(y) 


—  OO 


ty  <  OO 


The  integrals  in  (2.93),  (2.94)  can  be  computed  once  the  functions  K^T),  K2(T),  kJf(T ), 
"(T)  and  the  boundary  values  «*0,  a1  have  been  expressly  prescribed. 

We  can  also  invert  the  equation  (2.76)  in  order  to  find  explicitly  Ts(t)  =  $ _  1(t), 
0  <  t  <  tf.  where 
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Ts(t ) 

(2.95)  *(Te(<))  =  [  K°(1  ~  us(y))kfj{y)e 

T's(0)  M2K2(y)-ckff(y) 

Let  us  now  examinate  the  second  region. 

U)b  (tt0)al)6/3r 

We  first  discuss,  as  in  the  previous  case,  the  solvability  of  the  equation  (2.89)  for  the 

.mtral  value  Ts( 0).  It  must  hold,  as  in  case  ii)a,  condition  (2.86),  otherwise  2f.(o)  would 
be  negative. 

On  the  other  hand,  according  to  (2.67)  and  (2.69)  it  must  be 

(2.96)  Ts(0)<t£(c,). 

From  (2.95)  and  (2.96)  we  deduce 

(2.97)  Tbs(a0)  <  T*s(c,y, 

this  condition  is  equivalent  to  (cfr.  (2.48)) 

(2.98)  aQ  >  a0  =  ~YjJ  kff(n)dri 

0 

Examinate  now  the  equation  (2.89).  Formulas  (4.2.90)  are  replaced  by: 

(2.99)  B(Ts(a 0))  >  b,  B(Tg(c/))  >  0,  B'(s)  <  0  for  s  <  Tg(ct) 
where  B(s )  is  defined  by  (2.89). 

Therefore,  there  exists  a  unique  solution  s  =  Ts( 0)  of  (2.88)  if  and  only  if  B(Tl(c,))<b 
But  hy  the  definition  of  r|(>0)  (2.17)  we  have 


TUc,) 


4(«o) 


s<n(c,))"  "i^/ 

0  0 

Thus  provided  that  (2.97),  there  exists  exactly  one  value  Ts( 0)  such  that  (2.82)  is 
satisfied.  ’ 
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Whenever  (a^eq) e and  a0  >  a0  the  initial  situation  is  consistent  with  (2.66)-(2.68). 
Indeed,  we  have  25(0)  <  0,  since 


(2.100)  Tg(c)  >  T*s(c,)  for  c,  <  c  <  K2(0)/ku. 


If,  on  the  contrary,  c>K2(0)/ku  (that  is  there  is  not  any  T*s(c)  solution  of  (2.15)),  it  still 
holds  z5(0)<0.  Moreover,  we  have  zF(0)  >  0,  and  equality  holds  if  and  only  if 
rs(0)  =  Tj(c,). 

From  (2.74)  and  (2.75)  we  see  that  Ts(t)  >  0  whenever  T5(0)  <  Ts(t)  <  T*s(c,).  The  front 
speed  zs(t)  does  not  vanid=sh  as  long  as  Ts(t )  <  T*s(c,)  ( zs(t )  would  be  zero  if  and  only  if 
-FsM =  Ts(c)>  but  Pi  (2.100)  holds).  The  solution  Tg(t)  of  (2.74)  can  not  have 
asymptotes  t  <  T*s(c):  indeed,  the  quantity  on  the  right-hand  side  in  (2.74)  would  not 
vanish  if  T s(t)  tended  to  £.  On  the  other  hand,  the  derivative  T5(0  is  bounded  by: 


|is(<)|< 


(^2(Q)  +  O(^o)2 

^(1-KTS(0))) 


se  rs(o)<r5(0<n(C|). 


Thus,  the  solution  Ts(t)  of  (2.74)  must  reach  the  value  T$(c ,)  in  a  finite  time  7.  At  that 
time  we  have  zF(t)  =  0  and  the  unfrozen  region  of  the  soil  becomes  exhausted. 

Once  the  temperature  Ts(t)  has  been  calculated,  by  inverting  (2.76),  we  can  achieve  the 
boundaries  zF(t )  and  z:s(t)  by  means  of  (2.69)  and  (4.2.70),  respectively.  Since 
Ts(t)  <  Tg(c[)  for  0  <t  <7,  we  see  that  the  constraints  (2.66)  respected. 

By  (2.61)  we  find  the  final  thickness  of  the  frozen  fringe: 

Ts(i) 

Zs(I)=~IW0J  >  ° 

0 

Giving  a  geometrical  interpretation  on  the  (a^aO-plane,  as  we  did  for  the  set  <$0,  we  can 
say  that  the  left-side  boundary  of  the  region  Tj  is  the  vertical  straight  line 

noi) 

(2.101)  ao=-I-±^J  kJf(v)dn, 

0 

while  on  the  left-hand  side  'Jj  is  bounded  by  the  straight  line  rc  . 

When  t  =  0  we  get  the  condition  (2.99).  As  t  increases,  the  straight  line  moves  on  the 
right  towards  the  point  an^  when 

kfWd*> 
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the  isotherm  T  =  o  reaches  the  base  of  the  soil. 


n  he  ground  of  the  analysis  of  cases  i,>  and  i,>  introduced  in  proposition  2.2,  we  see 
that  tt  ,s  more  appropriate  denote  the  regions  and  5,  by  %(,,,  and  by  5lfe(i)). 

e  con  rary,  uring  the  process  of  lens  formation,  the  boundaries  of  III,  keep  at 

rest  Tins  corresponds  to  the  fact  that  the  solution  for  lens  formation  is  global,  while  the 
solution  for  frost  penetration  is  local. 

In  order  to  understand  better  the  deference  between  the  teo  cases  ,,  =  cc  when 

°  °'  prOpOsiti011  2  2>’  ™  I*™  to  make  some  remarks  about  the 

curves  C,  and  C2  defined  by  (2.37)  and  (2.38),  respectively. 

Consider  the  two  famil.es  of  curves  C,(»)  and  C#,  depending  on  the  parameter  6  =  r5(0) 


Pl(b)  rcr\C1,0<c<cl  P%(b)  =  rcriC2,  0  <  c  <  c, 

Where  the  the  straight  lines  rc  are  defined  by  (2.23).  The  points  p((t)  and  p|(»)  m  weu 
efined  according  to  the  properties  of  the  curves  c,  and  e2  (see  remark  2.1). 
h  <  b2,  it  is  easily  checked,  by  examinating  (2.37)  and  (2.38),  that 

(2.102)  |  P{(b2)  -o\<  j  Pffa)  -of  <  |  Pl(b2)  -Oj<  Ipc(6i)  _  0|,  0<c<c, 

and  equality  holds  if  and  only  if  c  =  C/. 

We  remark  that  the  slope  of  the  straight  fines  r0  and  which  fix  the  boundaries  of  the 

Zt  ccTr  a?  °U  l'  Hen“’  iMqUali‘ieS  (2102)  *>“*  -  •  -eases  the 

urves  e,W  and  C2(t)  move  towards  the  origin,  but  keeping  inside  the  angle  A. 

The  property  (2.102)  is  showed  in  fig.  4.2.  C 
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figura  4.2:  the  curves  C,(6j)  and  C,.(62),  i  =  1,2,  with  b1  <  b2. 


Let  us 
time  t  ( 

(2.103) 


take  now  and  the  corresponding  solution  of  (Sf,)  for  t  e  [0,<y),-  for  each 

=  [0,< y)  consider  the  curves 


2(2s(0)  =  {(o0,Qi):o-  +  c^^fzs(t)  + 


1 

kua0 


T’s(c) 

I  ^ 


2 W-ckfjW 


K,{V) 


dr)  =  0,0  <  c 
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According  to  (2.102),  the  distance  between  a  point  («„,«,)  fixed  in  %  and  the  curve 
V*s(<))  decreases  as  t  increases. 

Moreover,  zs(t)  =  0  if  and  only  if  (a0,aj)  e  C2(zs(*));  in  that  case,  Ts(t)  =  r*(c). 

Now,  if  t}  =  oo,  then  zs(t)  never  vanishes  and 

liZ o  M*VsW,t5(o)  =  (4K,«X),  4K>ai),n(c(a0, «!))). 

This  is  equivalent  to  state  that  the  curves  C2(^(0)  never  cross  the  point  (tt  *  ) 
corresponding  to  the  assigned  boundary  data. 

If,  on  the  contrary,  tf  <  oo,  we  have  (<*„,*,)  €  C 2(zs(tf))  and  zs(tf)  =  0 

Consider  the  solution  Ts(t)  of  (2.74)  for  «  >  V  We  have  Ts(lf)  „  0,  ,f  t<(<)  >  0  in  some 

ng  interval  ,  >  «/t  we  get  a  contradition  with  (2.66),  since  Ts(,)  would  be  greater  than 

sM  Ts(tf)  (thus  ckIt(Ts(t))  >  Ayjyijj),  while  .e(Ts(,n  defined  by  (2.75)  is  negative. 

n  the  other  hand,  if  Ts(t)  <  0  in  some  right  interval  of  tf,  the  solution  of  (2  74)  is  not 
consistent  with  (2.66),  (2.73).  1  j 

We  conclude  that  it  is  no.  possible  to  get  a  solution  of  (Sf,)  consistent  with  the 
prescribed  conditions  when  t  >  t ^ 

’  Pr0CeSS  °f  Ie°S  f0rma,ion  8tart  <“  <  =  «,•  As  a  matter  of  fact, 
s  )  7  S'  h  ‘lJtl  =  4(«o,o,)  is  solution  of  (Sji)  for  < >  (,  with  t  =  y  =  , 

The  time  <  =  is  a  transition  time  from  the  process  of  forst  penetration  to  the  other  of 
lens  formation.  We  remark  that 

&«&i>  „  <>, 

dz  — 

tWgtonX  P°int  <<’0’<,l)  ®iVe°  ^  ‘he  PreSlibed  da‘a  ‘ieS  OI  ‘he  boun<iary  e^s(‘,»  of 

Summarizing  the  results  got  in  parr.  2.1  and  2.2,  e  may  conclude  that  for  any 
assignment  of  the  data  («„,»,)  on  the  plane  and  of  the  initial  value  %(o)  =  b  <  a,  we  are 
in  position  to  foresee  if  a  process  of  lens  formation,  or  frost  penetration  or  none  of  them 

”  t'/r  “  SUffiden‘  Cl,eCk  if  the  P™4  K,«i)  lies  in  m  or  in 
T(fc)  -  T0(&)  U  %(b)  or  in  none  of  these  regions. 

The  following  proposition  sums  up  the  developed  analysis. 
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PROPOSITION  2.3. 


0 

ii)  If  3  T a  <  0  such  that 


Vj?  6  ( —  oo,0);  then  i.(d>)  =  =  0. 


<?  + 


T <r 

[  ML 

J  «i(v) 


dr)  -  0, 


then  l  and  <5  are  two  contiguous  regions  on  the  (o^-plane,  having  the  curve  C2, 
defined  by  (2.38);  as  a  boundary  in  common. 

Moreover,  whenever  (<*„,«,)  el(b),  the  system  (S^)  has  a  unique  solution  for  any  time 
t  >  0  such  that  the  temperature  Ts  is  constant,  as  well  as  the  thickness  of  the  frozen 
fringe,  the  growth  of  the  upper  boundary  of  the  soil  is  linear  with  respect  to  time.  When 
(case  «>  of  proposition  2.2),  the  system  (s£)  has  a  unique  solution  that 


attains  m  a  finite  time  (in  that  case  we  have  a  transition  to  the  process  of  lens 
formation)  or  infinite  (in  that  case  the  solution  is  global)  the  stationary  values 
introduced  in  the  proof  of  thr  proposition  2.2. 

When  (a0,ai)€V1  (case  ii)b  of  proposition  2.2),  the  system  (S ft)  has  a  unique  solution 

whose  mam  feature  is  that  the  isotherm  zF{t)  =  0  reaches  in  a  finite  time  the  base  of  the 
soil. 


The  regions  1(b)  and  «F(6)  are  outlined  in  fig.  4.3. 
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fig.4.3:  the  regions  £(6)  (lens  formation)  and  5(6)  (frost  penetration)  on  the  (a0?a,)- 
plane.  The  curve  C2(6)  is  a  common  boundary  to  the  two  regions;  the  region  5(b)  is 
bounded  on  the  left-hand  side  bythe  straight  line  a0  =  50(fc).  The  two  parts  50(6)  and 
^(6)  are  related  respectively  to  the  cases  ii)a  and  u)b  of  proposition  2.2. 
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3  Boundary  thermal  fluxes  depending  on  time 


We  are  going  now  to  discuss  the  solvability  of  (S/z)  (defined  in  par.  1.5),  together  with 
the  conditions  (C)  (par.  1.7)  and  the  constraints  (ffl)  (par.  1.8)  removing  the  assumption 
a0,  constant. 

Our  main  purpose  is  investigating  the  possibility  to  pass  more  than  once  from  one 
process  to  the  other,  in  order  to  simulate  a  process  of  penetration  of  the  front  zs,  which 
occasionally,  under  appropriate  conditions,  stops  giving  rise  to  the  process  of  lens 
formation.  In  the  previous  section  we  have  already  dealt  with  one  case  of  transition 
process  form  frost  penetration  to  lens  formation:  more  precisely,  it  occurred  when 
(<*o>ai)  6  %  and  tf<oo.  We  find  it  convenient  to  write  the  set  of  equations  (S/z)  +  (C)  +  (ffl) 
in  the  following  way: 

TS(t) 

\  zF(t)  -zs(t)  =  1-L^J  kff(r))dr) 


(S //)  +  (€)  + (A) 


?«>(*)  —  (i  -  + 


kfal(t)  -*uQo(<) 

Lpw 


Pih(t)  =  Pw< !«,(*)  +  e(l  -  vs){Pi  -  pjzs{t) 
Ts(i) 

zs(0)  =  6,  zT( 0)  =  H  >  b 

- Q~z -  *s(*)  =  0 


*s(<)<0  , 

dz 


>  0 


0  ^  z ^  z <  Zrp(i^ 


?  JO  >  ° 
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Moreover,  we  know  that  when  %(<)<0  (henc  ^»(g(l)’l)  =  0),  then  (s„)  +  (€)  +  , 

«cnbes  the  process  of  frost  penetration;  when  is(()  =  0,  the  process  is  lens  formation 
We  read  the  assigned  data  a0(t),  ax(t)  as  a  curve 


(3-1)  7  = 

which  is  prescribed  on  the  («0,  ..(-plane,  with  „J«)  >  0,  .?(i)  >  0,  ,  6  (0ig. 

Usmg  the  same  notations  as  in  sect.  2,  we  will  denote  by  lW  and  5(»)  the  regions  on 
e  plane  where  the  initial  conditions  allow  the  occurrence  of  lens  formation  and  frost 
penetration,  respectively,  according  to  the  summary  given  in  proposition  2.3.  We  will 
assume  that  the  temperature  T,  verifying  (2.29)  exists,  so  that  it  is  1(4)  ^  „,  5(J)  1 

3.1  Preliminary  results 

1)  If  7  C  1(b)  Vi  e  [0,to],  then  (Sfl)  +  (C)  +  (A)  has  a  unique  solution  describing  a  process  of 
ens  formation,  such  that  the  temperature  Ts,  the  thickness  of  the  frozen  fringe  b 
ana  the  hydraulic  flux  qw  depend  on  time . 

Actually,  call 


C7  (a?(i),»7(i))  = 


LpwK*m 


or,  more  briefly,  c7(i)  and  assume,  contraryto  our  claim,  that  a  solution  of  (S„)  +  (C)  +  (fl) 
venfies  zs(r)  <  o  at  some  time  r  6  [0,io],  By  virtue  of  (2.72),  we  would  have 

(3-3)  K2(ts{t))  <  c7  (r)kff(Ts(T)). 

We  denote  by  TJM  the  temperature  satisfying  hT2(TJ( r))  =  ,-|,W))  (il  is 
umvocally  determ, ned  since  WM,.?<r()  e  A  thus  <  *,(„)/»„  and  (2.15  has 

exactly  one  solution).  From  (3.3)  we  deduce: 

(3-4)  Ts(t)<T*s(t) 

Moreover,  from  (S/,)  +  (C)  +  (ffl)  we  easily  find: 
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(3.5)  6>,sW=  -JL  Z(Ts(r))>  -jJjU  Z(T's(r» 


where 


3(»)4'W| 


*!<») 


/ 

'♦/ 
V  o 


*//W 


*lO?) 


5 

rfj?  +  7^/  kff^drf 


(cfr.  (2.70)). 

Inequality  (3.5)  yelds  that  the  point  P(t)  =  (a^(r),a^(r))  lies  in  «F(*5(r))  and  that 
P(t)  £  C2(z5(r)).  So.  we  get  a  contradiction. 

We  conclude  that,  whenever  (atf(t), a?(t))  e  1(b),  it  is  possible  to  find  a  solution  of 
(sfl)  +  (C)  +  W  only  if  zs(t)  =  0,  t  e  [0,<o]. 

On  the  other  hand,  talcing  zs(t)  =  0,  t  e  [0,<o],  the  solution  of  (S/;)  +  (C)  +  (A)  is  for  each 
time  t 


(3.6)  ,,(t)=;WWiVgffl 

r5W 

(3.7)  zF(t)  =  ^  J  kjf(T])dri  +  6 

where  T5(t)  is  the  uniqur  solution  (cfr.  (2.13))  of 


(3.8) 


+  c^(t) 


*u«o(0| 


( 


Kn 


Ts(t ) 


b+kAm!  *"(VH+/ 


Ts(t) 


K^)-c^t)kif{r,) 

Klin) 


- dr}  =  0. 


The  speed  irr(t)  and  the  height  of  the  soil  are  given  by 

(3.9)  ir(<)  =  5?U*) 

.  _  .  /  fc,4(r)-£  ctf(r) 

(3.10)  zyW  —  ff  +  J  2^ - ^r 

0 

Finally,  we  achieve  the  temperature  T  and  the  pressure  pw  in  any  point  of  the  frozen 
fringe  and  at  any  time  t  by  means  of  (1.20)  and  (1.28). 

If  we  are  interested  in  checking  when  the  thickness  of  the  frozen  fringe  b  -  zF(t)  and  the 
temperature  T s(t)  are  increasing  or  decreasing  with  respect  to  time,  we  have  to  derivate 
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(3.7)  and  (3.8): 

/O  ,  /  ,7  Ts(t)  > 

(  ’  M<)  =  M  J 

v  n  ) 


rsW 

"  *35;  -  *(') /  ‘«w(  3^  -  jrLy)(, 

rUi\h  .  ( T  <4\\  TA  Znri  7TT\  T7 - — - - - - - 


(3.12)  Ts(t)  = _  n  V“° 

r^(t)kfJ(Ts(t))  K2(Ts(t ))  -  c7(<)it^(Ts(<)) 

A'o  W)i 

We 


IV 

make  the  following  remarks  (dropping  the  apex  y,  for  simplicity): 

■)  If  do(<)<0,  «i(i)>0,  then  fs(  t)  <  0,  i_(()  <  n  i  ms,  T„  ..  , 

(  oW)  +  >  °>  tie"  TM  <  «,  V(<)  <  o,  U<)  >  o  (curve  r,  in  fig.  4.4). 

(i)<0.  If  (a0(*))2  +  (aqft))2  >  0, 
UU11^  ocuac  ^curve  72  m  fig.  4.4). 

For  a  more  general  curve  7,  we  have  Ts(t)  <  0  if  and  only  if 

Ts(t) 

(3'I3)  ». 

Inequality  (3.13)  is  not  easy  to  solve  in  terms  of  since  Tj(<)  depends  on 

and  ajft)  by  means  of  the  implicit  relation  (3.8). 


U)  K  "oW-0)  ^<0,  then  Ts(t)>  o,  *F(*)>0,  *„(*)<  0.  Ii 
the  inequalities  hold  in  the  strict  sense  (curve  ?2  in  fig.  4.4). 
a  more  general  curve  7,  we  have  Ts(t)  <  0  if  and  onlv  if 
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fig.  4.4:  process  of  lens  formation  along  the  curves  7j  and  y2.  Along  7l  the  temperature 
at  the  base  of  the  lens  Ts  decreases,  the  thickness  of  the  frozen  fringe  and  the  hydraulic 
flux  are  increasing  with  respect  to  time.  Along  y2  T s  is  increasing,  while  the  thickness  of 
the  frozen  fringe  and  the  hydraulic  flux  decrease. 


2)  If  the  curve  7  is  such  that  (<))  €  7(z5(<))  for  any  time  t  €  [0,to],  where  <5( zs(t )) 

is  the  set  bounded  by  the  curve  C2(z5(t))  and  defined  by  (2.38),  then  (Sfl)  +  (C)  +  (A)  does 
not  admit  solutions  describing  lens  formation  ( that  is  with  zs(t)  =  o). 

In  order  to  check  2),  it  is  sufficient  to  argue  as  in  point  1),  by  assuming  ab  absurdo  that 
a  solution  (S//)  +  (C)  +  (A)  such  that  is(t)=0  may  exist.  We  easily  find  a  contradiction 
with  the  assumption  (a3'(t),o7(0)  e  5(z5(t)). 
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Let  us  keep  now  the  hypothesis  (aft*), «?(<))  e  «F(*5(t))  and  write  formulas  which 
generalize  (2.71)  and  (2.73),  in  order  to  encompass  the  present  case  of  thermal  fluxes 
depending  on  time: 


(3.14)  i 

s(t)-  (1-* 

(3.15)  is(t)=--A 

KuaC 

(  ( 

,  K0a\{t) 

[  ku(Tsm 

(l-^Uy/T^))' 


-c7(t) 


K0  d 


wtsmH  y\ 


H 


K2{Ts{t)) 


KpW\dTi\K2{Ts{t)) 


y,] 

a+  J  - dTl 


+ 


kff(Ts(t ))  L 


V  o 


Ki(v) 


Kn 


TS(t) 


s(t)  + 

\ 


+ 


Ko  J  k}S{r>)dr> 


) 


By  combining  (3.14)  and  (3.15),  we  get  the  following  ordinary  differential  equation: 

fi^(rs<(,)rs(«,  +  °y$HTsW))  =  <  W 

where  <p(s)  has  been  defined  by  (2.75)  and 


(3.16)  Eo(1~l'sK 
(k. 


(3.17)  V-(s)  =  ^4^ 


/ 


K2{s) 


ff+I 

\  o 


5  K2 


KM  dt> 


S 

~~ko  J  kfMdri- 


If  we  take  constant  a0,  a,  in  (3.15),  we  get  the  stationary  equation  (2.74). 
The  functions  <p(s)  and  ip(s)  are  surely  negative  fors  <  T^c,). 

Since  we  set  (<*#(<), a}^*))  e  ^(z^t)),  Vfe  have,  whenever  0  <  c  <  K2(o)/ku: 


(3.i8) 


/ 


\ 


+  J 


o 


KM 


-dr)  >  0 


On  the  other  hand,  if  Ts(t),  zs(t)  verify  (S/f)  +  (C)  +  (A),  it  must  hold  for  any  t  e  [o,<0]: 
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(3.19) 


<r  + 


K2(Ts(t))  ku*m 


( 


kff(Ts(t))  Ko 


TS(t)  \ 

5(<)  +  ^fe/  kfWdTtj  + 


,f(<)  0 

+  J  K^rj)  dT>-° 

0 


From  (3.18)  and  (3.19)  and  reminding  that  c7(t)  =  ,  we  deduce: 

t/ATs(c7m 

(3.20)  Ts(t)  <  T*s(c7(i)) 


and 

(3.21) 


kff(Ts(t)) 


If  c  >  iT2(0)/&u,  then  (3.21)  obviously  holds. 

We  may  conclude  that,  when  7  c  ?(^5(<))  (that  is  7  is  ever  on  the  right  side  of  the  curve 
C 2(zs(t)))i  f°rnI  zs(i)  ^as  property  zs(t)  <  0.  On  the  other  hand,  if  we  know  the 

solution  Ts(t)  of  (3.16),  we  can  achieve  from  (2.70)  (which  is  valid  even  in  the  case  a0, 
a1  depending  on  time)  the  boundary  zs(t)  and,  by  means  of  (2.103),  the  profiles  of  the 
curves  C2(z s(t))  that  delimitate  9:(z5(^))  at  each  time  t.  By  comparing  the  path  7 
assigned  as  in  (3.1)  with  the  regions  we  are  able  to  evaluate  how  long  the 

process  of  frost  penetration  lasts. 

If  we  add  the  hypothesis  &o(t)  <  0,  then  we  get  from  (3.16): 


(3.22)  Ts(t)  >  0. 


3.2  A  transition  process 


We  are  now  in  position  to  describe  a  process  where  the  curve  7  moves  form  the  region 
to  JL  and  viceversa. 

Consider  a  path  7  =  (aj^b^i  (Ob  t  e  [<M0],  such  that  (c$(0),a7(0))  e 
Assume,  at  least  for  the  moment,  that  the  additional  hypothesis 


(3.23)  e$(t)<0 
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holds. 

Getting  Ts  from  (3.16)  yelds 


(3.24) 


Ts(t)  = 


_  1  (  «7o  (t) 

«J(i) 


^s(<))  + 


-c^(t) 


Equation  (3.24)  shows  that 


[  /or  Ts(<)  <  min{r^(0),Tj(C/)},  0  <  c7  (t)  <  K2(0)/k 
(3.25)  fs(t)  >  0  | 

(  for  Ts(t)<T*s(c,),  c>K2(0)/ku 

Furthermore,  the  derivative  Ts(t)  is  bounded  by  a  constant  if 

sup  laJ(t)l<oo 

<e[o,<0]'  '' 

and  Ts(t)  can  not  have  an  asymptote  i  ^  T£(c7(i)). 

We  may  conclude  that,  whenever  y  c  ^i{zs(t))  for  any  t  e  [o,<0]  (curve  Tl  in  fig.  4.5),  then 
Ts(t)  reaches  the  temperature  T*s(c,)  in  a  finite  time  7;  in  that  case,  the  isotherm  T  =  0 
matches  the  base  of  the  soil,  that  is  zF(t)  =  0.  The  final  thickness  of  the  frozen  fringe  is 

TsO) 

(3'26)  ^  =  ku^0<J)  l  kff^dT}- 

0 

The  second  possibility  for  the  path  y  is  when  there  exists  a  time  tf  <t0  such  that 
n(tf)^(tf))eC2(zs(tf))  (fig.  4.6);  hence  Ts(tj)  =  T*s(c~<(tf))  and  zs(tf)  =  0.  We  remark 
that  m  the  present  case  of  thermal  fluxes  depending  on  time,  Ts(tf)  is  not  generally 
zero,  as  it  occurs  in  the  case  a0,  constant. 

It  may  happen,  as  we  saw,  that  Ts(t)  tends  to  T‘(c7)  in  infinite  time  (tf  =  oo).  In  such  as 
case,  aj  must  verify  =  0,  otherwise  it  is  impossible  to  have  stationary  solutions. 

Conversely,  assume  that  tf  <  oo  and  that  (a^WW)  e  in  some  interval  t>tf 

(curve  y2  m  fig.  4.6).  In  that  case,  a  lens  starts  growing  at  the  time  tf  and  at  the  height 

*  =  Zs(tf)'  As  long  as  7  lies  in  l(zs(lf));  which  has  fixed  boundaries,  the  thickness  of  the 
lens  goes  on  increasing;  the  temperature  Ts  and  the  thickness  of  the  frozen  fringe 
zs(i/)~zF(t)  change  with  respect  to  time,  according  to  the  statement  1)  in  par.  3.1. 
Essentially,  the  process  of  lens  formation  developes  following  two  directions: 
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i )  there  exists  a  time  ^  >  tf  such  that  (ao(*i)>ai  (*i))  e  ^(^(fy)):  in  that  case  the 
frozen  fringe  invades  the  whole  unfrozen  soil  and  zF{t^)  =  0  (branch  72  in  fig.  4.6); 


ii)  there  exists  a  time  t2  so  that  kuo$(t2)  —  k^aj(t2):  in  that  case,  the  water  flux  qw 
vanishes  at  t  =  t2  and  a  melting  process  will  take  place  for  t  >  t2,  if  7  goes  under  the 
straight  line  kua0  =  k^ar  (branch  j2  in  fig.  4.6). 

For  more  general  curves  7,  which  has  not  the  property  (3.23),  further  posiibilities  occur, 
in  addition  to  the  ones  just  introduced. 

First  of  all,  Ts(t)  may  also  be  negative,  if  0^(2)  >  0. 

As  long  as  Ts(t)  <  mm{T5(c7(*)),Tj(cj)},  the  velocity  of  the  front  z5(*)  and  the  boundary 
^  =  zF{t)  never  vanish.  If  Ts(t)  decreases,  it  tends  to  a  stationary  value  Tt  if  and  only  if 


(3.27) 


_  go(0 

H7(<))3 


2  K2 (Tf) 

uhf(Tiy 


Conversely,  if  Ts(t)  increases  up  to  the  temperature  Tj(c7(Z))  when  t  =  if  and 
7  C  JL(zs(tf))  for  t  >  tp  then  the  process  of  lens  formation  starts. 

In  addition  to  the  possibilities  stated  above  as  for  the  stop  of  the  process  of  lens 
formation,  a  third  case  may  occur,  when  the  curve  7  goes  into  the  region  ^(zs(tf))  once 
again,  say  when  t  —  t3  >  tj  and  a  second  process  of  frost  penetration  takes  place  with  the 
initial  condition  zs(tf)  for  the  boundary  zs(t)  (branch  72  in  fig.  4.6). 

Iterating  the  described  process,  one  can  get  a  sequence  oflenses  and  each  ice  layer  is 
separated  from  the  previous  one  (which  is  above)  by  a  layer  of  frozen  soil,  corresponding 
to  a  process  of  frost  penetration. 
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fig.  4.5:  development  of  the  process  of  frost  penetration  along  the  curve  Tj.  For  t 
isotherm  zF  reaches  the  base  of  the  soil. 


(ot'MJX o)) 


fig.  4.6:  process  of  frost  penetration  starting  from  the  boundary  fluxes  values 

(ag(0),^(0))€5(2s(0)). 

For  t  =  tj  the  velocity  of  the  front  zs  vanishes  and  lens  formation  succeeds  frost 
penetration.  Along  the  branch  the  process  stops  since  Zp  reaches  (by  decreasing)  the 
base  of  the  soil  when  t~t^  along  7*  the  water  flux  qw  vanishes  when  t  =  t2 ■  along  7^  a 
second  process  of  frost  penetration  takes  place  when  t  =  t3. 
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4  Temperatures  specified  at  the  boundaries 


We  write  again  the  set  of  equations  (Stmp)  +  (C)  +  (A),  defined  in  sect.  2  that  describes 
t  e  freezing  process  in  the  case  that  we  assign  the  temperatures  at  the  extremities  of  the 
soil  h(t)  =  T(0,t)  and  g(i)  =  T(zT(t),t),  h(t )  >  0,  g(t)  <  0: 


Ts(t) 


(Simp)  +  (C)  +  Ws 


,  (4-i)  / 


.mo 


(4.2)  Lfw,M  =  (1  - „s)ePwL i  (,) _ i  _  W) 

/arrW-*sW  *F(0 

(4.3)  PtzT(t)  =  pwqjt )  +  e(l  -  us)(Pi  _  Ptu)i5(<) 


(4.4)  +  / 


T <j(t)  if  (  ^  Qt /  \ 

,  r  R2{T,)~~rmrkffM 


*M 


-dr) 


_  ffJO 


*f(0 


(4-5)  +g,(^,W>s(()) 

(4.6)  *5(0)  =  6,  zT( 0)  =  H  >  6 


ta  71)  5P«,(*5(<)»<)  ....  _ 

- Qi - ^(<)  =  0 


(4.8)  is(t)  <  o  , 


(4.9) 


dpJSsJ) 


dz 


-  >0 


(4.10)  0  <  zF(t)  <  zs{i)  <  zT(t ) 

C4-11)  ?J0>  0 


The  unknown  quantities  ate  the  water  flux  ,„(<),  the  boundaries  ,#),  %(()  and 
the  freezing  temperature  TM  and  the  water  pressure  p.fed).,,.  We  recall  again  that 
w  en  zs(t)<  0  (in  ttha  case  the  water  pressure  gradient  vanishes  at  z  =  zs,  see  ...),  we 
have  a  frost  penetration  process;  when  zs(t)  =  0,  lens  formation  occurs. 


Remark  4.1.  In  theory,  we  should  make  use  of  the  results  in  sections  2  and  3,  in  order  to 
solve  the  present  case,  basing  ourselves  on  the  relations 


a0(t)  = 


«i  (t) 


zT(t)  —  zs(t) 


Considering,  for  istance,  the  case  of  lens  formation,  once  h{t)  and  g(t)  are  given,  we 
should  find  aQ  and  a1  such  that  the  following  equations,  coming  from  (4.1),  (3.8)  e 
(3-10),  are  satisfied: 


a0(t)  = 


Ts(t) 

KKi)-j  kff(’n)dv 

o 


aa(<)=- 


Ts(t)-g(t) 


H  +  TiJ  (kfa  i(r) _  kuao(T)YT 


cr  +  c(t) 


l'uo-o(0| 


A'n 


b  + 


Ts(t ) 


+ 


Ts(t ) 

/ 


A2(t?)-c(f)^//(7?) 


t/77  —  0 


where  c(t)  is  defined  by  (3.2). 

The  system  we  have  just  written  is  anything  but  simple,  mainly  owing  to  the  implicit 
dependence  of  Ts  on  the  fluxes  a0  and  ax  in  the  third  equation. 

Therefore,  we  believe  it  is  more  convenient  and  more  interesting  to  solve  directly  the 
problem  (4.1)- (4.11). 

In  paragraphes  4.1  and  4.2  we  will  discuss  separately  the  cases  of  lens  formation  and 
frost  penetration,  respectively.  Finally,  we  will  investigate  on  the  posssibility  of  a 
transition  process  (par.  4.3). 

4.1  Lens  formation 

From  (4.1)  we  get  zF(t)  in  terms  of  Ts(t): 
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(4.18)  ¥>i(s,*)  >  0,  ^(s,*)  <  0  per  s  <  0,  t  >  0. 

From  (4.13)  and  (4.14)  one  gets  the  following  identity  for  zT(t)  as  a  function  of  the 
temperature  Ts: 

t4  1Qi  _ i. 

T  ~  +  LP,  Vii't'sW)  +  ViC11  s(t),i)<P3(Ts(t),t)  foW ~  Ts(*))  =  <P(Ts(t),t). 

Denvating  (4.19)  with  respect  to  time  and  comparing  with  (4.13),  we  get  the  following 
ordinary  differential  equation  for  Ts(t ): 

(4.20)  _l2(Ts{t))  _  d<p(Ts(t),t)  d<p(Ts(t),t) 

1  '  v>i  (Tsm~  dTs  T^)+ — oi — • 
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We  must  add  an  initial  condition  for  T5,  that  is  not  esplicitly  prescribed  by  the  set  of 
equations  (4.1)-(4.11):  T5(0)  has  to  be  computed,  using  the  known  initial  values  ft(0), 
#(0),  H  e  6,  and  solving  the  equation  y?(T5(0),0)  =  H ,  that  is 


(4.21) 


jfaggW)  ■ 

^,(t5(0),0)  - 


kf  g(0)  -  Ts(0) 
LPi  H-b 


+  ^3(Ts(  0),0) 


Our  first  purpose  is  to  discuss  the  solvability  of  (4.21).  To  tis  end,  we  rewrite  the 
conditions  (4.7)-(4.11),  which  in  case  of  lens  formation  become 


(4.22)  >0 

(4.23)  0  <  zF(t)  <  b 


(4.24)  zT(t )  >  0, 


in  terms  of  the  temperature  Ts .  We  will  first  concentrate  our  attention  on  the  initial 
time  /  =  0:  we  will  look  for  suitable  conditions  on  the  initial  data  h( 0),  ^(0),  H  and  6,  so 
that  the  initial  situation  is  consistent  with  the  conditions  (4.22)-(4.24),  computed  for 

We  start  with  the  following 

LEMMA  4.1.  If  there  exists  a  temperature  T ^  such  that  <p2(T^)  =  0,  then  for  all  t>  0 
there  is  exactly  one  value  Tp(h(t))  satisfying  the  equation 

(4.25)  <P2(Tp(h(t))  =  ^(Tp(A(0)^)^(rp(ft(<)),<) 

Dim.  We  write  (4.25)  as  follows: 


(4.26)  kuh(t )  =  r(Tp(A(t))) 


where 

(4.27) 


r(*)  = 


( 


.  kff(s) 

0*2(*) 


r 


S  K2(v) 

0 


KM 


dy 


\ 

/ 
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Let  us  consider  the  function 


*■(«)  =  '+  J 


■s 


KM 


dr},  s<  0. 


From  the  hypothesis  (Hu)  and  from  (1.31)  one  finds: 

(4.28)  F'(s)  =  J l^Mdr)  >  F(g)  >  VM  for  5  <  0,  F( 0)  =  ?>2(0)  =  <r. 


Assume  ab  absurdo  that 


(4.29)  F(s)  >0  Vs  <  0. 


By  virtue  of  the  properties  of  the  functions  K2  e  kjf  we  have  just  recall,  for  any  fixed 
s  <  0  there  exists  exaclty  one  value  c  e  (0,/sT2(0)/y  such  that  (see  also  (2.15)) 

&2(s)  ~  ckff(s)' 


If  (4.29)  is  valid,  it  would  follow 


Tj(c) 

<4'3°)  '  + 1 


but  (4.30)  is  in  contradiction  with  Lemma  2.2  of  sect.  2.1. 
Therefore,  there  exists  s,  <  0  such  that  F(Sl)  <  0. 
Reminding  (4.28)  and  (1.31),  we  find,  for  s  <  Sl: 


kff(s) 


kff(s) 


«  =  ~  K° >  ~  K°1^FM  >  “ 


°K2(s) 


°K2(s ) 


Since  K2(s)  tends  to  zero  for  s  going  to  -  00,  we  conclude  that  the  function  r(s)  reaches 
somewhere  the  value  kuh(t ):  in  other  words,  (4.26)  has  at  least  one  solution  for  each 
<  >  0  (notice  that  r  does  not  depend  explicitly  on  time  t). 

The  uniqueness  of  the  solution,  that  we  call  T„(A(t))  to  point  out  the  dependence  on 
time  only  through  the  boundary  temperature  h,  is  immediately  achieved  by  observing 
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that  r(0)  <  0  and 


>0  for  Ta  <  s  <  0 

(4.31)  r'(s)  <  —  0  for  s  — 

<0  for  s  <  T a 

Finally,  taking  into  account  that  r(s)  achieves  its  negative  minimum  for  s  =  T we  have 

(4.32)  Tp(h(t))  <  Ta.  □ 

The  following  result  allows  us  to  estabilish  when  the  equation  (4.21)  has  soluctions 
consistent  with  the  conditions  (4.22)-(4.24),  computed  for  t  =  0. 

PROPOSITION  4.1. 

i) If 

(4.33)  j  +  o'  >  0  Vs  <  0, 

o  1 

then  (4.13)-(4.14)  has  not  solutions  consistent  with  (4.22)-(4.24). 

ii)  Suppose  that  there  exists  one  temperature  T a  such  that 

T* 

(4'34)  / 

0 

then  the  following  cases  occur: 

ii)a  if 

(4.35)  kfg{ 0)  >  kfT„  +  (H- b)LPi<p3(T„, 0), 

then  (4.21)  has  not  solutions  consistent  with  (4.22)-(4.24)  computed  for  t  =  0. 
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«>  If 


(4.36)  kfTp  +  (tf  6)Xp^3(Tp,0)^l  +  <  kfg{ 0)  <  kfTa  +  (ff- bfip^T,, 0) 

theTl  the  initid  data  are  consi^nt  with  (4.22)-(4.24)  <md  (4.21)  »  aata/ierf  /or  exactly 
one  value  Ts( 0)  e  {Tp(h(0),g(0),b,H),  Tj  where 

(4-37)  Tp  =  max{Tp(h(0)),Tm(h(0),g(0),b,H)}, 


while  the  temperature  Tp(h( 0))  »  defined  by  (4.26)  {computed  for  t  =  0) 
one  value  such  that 


and  Tm  is  the 


(4.38) 
ii)c  If 


kM*)-TJ_ 
LPi<P3(Tm, W~H~b 


(4.39)  kfg(0)  <  kfTp  +  (H  -  b)LPi^(Tp^l  +  PwL^^jy 
then  we  conclude  as  in  ii)a. 

D,m.  From  (4.1)  and  (4.23)  it  follows  that  the  temperature  Ts  cannot  be  positive. 
Hence,  if  Ts(t)  is  the  solution  of  our  problem,  we  get  from  (4.18): 


(4.40)  Vl (Ts(t),t)  >  0,  v3(Ts(t),t)  <  0 


Assume  that  the  hypothesis  (4.33)  holds:  this  means  that  <p2(s)  >  0  for  all  s  <  0  We  see 
immediately  that  the  equation  (4.13)  cannot  be  satisfied  by  a  function  *,.(,)  such  that 
condition  (4.24)  is  fulfilled.  So,  the  case  i)  is  proved. 

Let  us  now  assume  that  there  exists  a  temperature  T„  which  depends  only  on  the 

r2  and  °n  ^Ut  DOt  °n  the  b0Undary  temperatures  h  and  g)  such  that 
(4.34)  holds,  that  is  <p2(T„)  =  0. 

Since  <p2(s )  is  increasing  for  s  <  0,  we  have: 


(4.41)  <p2(s)  <  0  if  and  only  if  s  <  Tg 
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Therefore,  the  solution  Ts(t)  has  to  verify: 


(4.42)  Ts(t)  <  T„ 


otherwise  (4.13)  is  not  consistent  with  (4.24). 

On  the  other  hand,  from  (4.21)  and  from  (4.40)  deduce  that  it  must  hold 


(4.43) 


zt(0)  =  H  <b  + 


kfim-Tsm 


Fi(Tsm  ■ 


We  remark  that  Fj(s)  is  an  increasing  function  for  s  <  0.  So,  if 


(4.44) 


H>F  i(0)  = 


kfg(0)\ 

Km)' 


no  initial  value  T5(0)  will  satisfy  the  condition  (4.43). 
On  the  contrary,  let  be 


(4.45) 


H  <b[  1 


kfg{ o)\ 

Km)] 


we  easily  see  that  there  exists  exctly  one  value  Tm,  g( 0)  <  Tm  <  0,  depending  on  h( 0)  and 
on  y(0),  such  that  =  H.  We  notice  that 


(4.46)  Tm  >  ff(0) 


(indeed  F^g)  =  b  <  H  =  FX{TJ). 

Owing  to  the  properties  of  F11  we  can  conclude  that  (4.43)  is  verified  if  and  only  if 


(4.47)  Ts( 0)  >  Tm. 


By  comparing  (4.42)  computed  for  t  —  0  with  (4.47)  we  see  that  it  must  hold: 


(4.48)  Tm  <  T,. 

Condition  (4.48)  is  equivalent  to  F(T c)  >  H,  that  is 
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(4.49)  kf{g{ 0)  -  Ta)  <{H  -  b)LPi<pz{Ta, 0). 

Therefore,  if  (4.35)  holds  we  see  that  no  initial  value  T5(0)  is  such  that  (4.43)  is 
satisfied.  Thus,  the  case  ii)a  of  proposition  4.1  is  proved. 

Assume  now  that  (4.49)  holds. 

We  have  still  to  impose  the  condiotion  on  the  water  pressure  gradient  (4.22),  that  in 
terms  of  the  temperature  Ts  becomes 

(4-50)  <p2(Ts(t))  >  ^))  ^{Ts{t),t)^{Ts{t),t) 

Condition  (4.50)  is  equivalent  to  (see  (4.27)) 

(4.51)  r(Ts(t))  <  kuh(t). 


By  virtue  of  Lemma  4.1,  whenever  (4.34)  is  true  there  exists  exactly  one  temperature 
Tp(h(t))  such  that  (4.51)  vanishes. 

Recalling  (4.31),  we  see  that  (4.51)  is  true  for  a  time  t  >  0  if  and  only  if 

(4.52)  Ts(t )  >  Tp(h(t)). 


Now,  the  two  posiibilites  can  occur  (see  (4.32)): 

2)  ^(^(O),!?^))  <  Tp(h( 0))  <  Ta 

Let  us  write  again  (4.21)  in  the  following  way: 

(4.53)  4>(rs(0))  =  F1{TS{  0))  -  H 

where 


(4.54) 


*(s)  =  j,(s,0)  =  (H-b) 


fg(f) 

AiM)y3(S,0)  > 


s  <  0 
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and  F^s)  is  defined  by  (4.43).  We  have: 


(4.55) 


(4.56) 


Thus: 


*(*)  >  0  for  s<Ta,  9(T„)  =  0 


(4.57)  tf'(s)  <  0  if  Tp  <  s  <  0,  *'(Tp)  =  0. 

In  the  case  1,  the  equation  \£(s)  =  F1(s)-H  has  exactly  one  solution  s  =  Ts( 0),  without 
further  conditions,  since  Fx(Tff)  >  H  (see  (4.49)),  F1(Tm)  =  H,  Fj(s)  >  0  for  (at  least) 
Tm  <  s  <  0.  Moreover,  we  have  Tm  <  Ts( 0)  <  Ta  and  Ts{ 0)  =  Ta  if  and  only  if  Tm  =  Ta. 

We  point  out  that,  by  virtue  of  the  properties  of  Fx(s)^  one  has 

(4.58)  Tp  <  Tm  if  and  only  if  kfg( 0)  >  kfTp  +  (H-  b)Lp^pz(Tp} 0). 

We  conclude  that  in  the  case  1)  we  find  exactly  one  initial  temperature  T5(0)  consistent 
with  the  prescribed  constraints  if  and  only  if 

(4.59)  kfTp  +  (H  -  b)LPi<p3(Tpfi)  <  kfg( 0)  <kfTa  +  {H-  b)Lp#z{T„, 0). 

The  inequality  (4.59)  is  well  defined  by  virtue  of  (4.52)  and  of  the  increasing  profile  of 
the  function  y?3(s,* )  with  respect  to  s. 

In  particular: 

kfg( 0)  =  kfTa  +  (H-  b)Lpiip3(TtT,Q)  if  and  only  if  T„  =  Tm  =  Ts( 0); 
kfg(0)  =  kfTp  +  (H-b)LPi<p3(Tp,0)kfg(0)  if  and  only  if  T p  =  T m. 


In  the  case  2),  the  equation  *(«)  =  F(s)-H  has  exactly  one  solution  if  and  only  if 
(4.60)  nTp)>Fx(Tp)-H. 

In  that  case,  the  solution  is  in  the  interval  [Tp,T„)  and  T5(0)  =  Tp  if  and  only  if  in  (4.60) 
is  true  the  equality. 

Taking  into  account  that 

(4-61)  *(: Tp)  =  (H-b)PwL 

k/A1  p> 

we  see  that  (4.60)  is  equivalent  to 

(4.62)  t/S(0)  >  k,Tr  +  (H-  l)iW3(y|  (i  +  PwL^^j  . 

Furthermore,  kepping  in  mind  (4.58)  once  again,  we  see  the  in  the  case  2)  we  find 

exactly  one  initial  value  for  the  freezing  temperature  Ts  consistent  with  the  imposed 
constraints  if  and  only  if 

(4.63)  kfTp  +  (ff-  b)LPi<p3(Tp,0)  ^1  +  <  kfg(Q)  <  kjTp  +  {H_  b)LPitp3{Tpfi). 

Condition  (4.63)  is  well  defined  by  virtue  of  (4.18). 

In  particular,  we  find: 

(4.64)  kfTp  +  (H  —  b)LPi<p3(Tp, 0)  ^1  +  PwL^^j  =  kfg{ 0)  se  e  solo  se  Ts( 0)  =  Tp. 

Putting  together  (4.59)  and  (4.63)  and  e  defining  Tp  as  in  (4.37),  we  achieve  (4.35)  and 
we  conclude  the  case  ii)b. 

We  remark  that,  by  the  definition  of  Tp  and  by  (4.46),  we  have 

(4.65)  Ts(0)  >  g( 0). 

From  (4.60)  we  deduce  also  that  if  (4.39)  holds,  no  initial  value  Ts(0)  is  suitable  in  order 
to  have  (4.12)  (computed  for  !  =  0)  satisfied.  Thus,  also  the  case  „>  is  proved. 
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Remark  4.2.  Since  dr{Tp(h))/dh  >  0,  we  see  that  if  h(t)  is  increasing  (decreasing),  the 
temperature  Tp  decreases  (increases). 

Our  next  aim  is  to  disuss  the  solvability  of  equation  (4.20). 

We  will  assume  from  now  on  that  condition  (4.36)  is  satisfied,  otherwise  the  initial  data 
do  not  allow  to  get  solutions  consistent  with  the  prescribed  constraints. 

We  consider  at  first  the  simper  case  in  which  the  boundary  temperatures  are  constant: 

(4.66)  h(t)  =  h>  0,  g(t)  =  g<  0. 

We  remark  that,  in  that  case,  the  functions  <px  and  <p3  do  not  depend  explicitily  on  time 
and  the  equation  (4.20)  reduces  to 

(4-67)  ~Sra=v'(Ts(,))}>(1) 

(for  simplicity  of  notation,  we  omit  the  second  argument  of  the  functions  <p1,  tp3  e  tp] 
moreover,  the  apex  for  tp  denotes  the  derivative  with  respect  to  Ts). 

Integrating  (4.67),  one  gets 

Ts(t) 

(4-68> 

*s(0) 


We  compute  explicitily  tp'(s)  making  use  of  (4.19): 


(4.69) 


#>'(*)  =  ~ 


kJ_  H^b 

LPi  ^O) Ws)  +  B-b) 


*'(*) 


where  ^  is  defined  by  (4.54). 

From  the  definition  of  tp3,  we  have: 


y>3(s)  <  0,  <p3(s)  - 


kff(s) 

bLPi 


>0 


for  s  <  0 


Furthermore: 
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i>(s)  +  H-b>0  per  s  <  T0,  i>'(s)<  0  forTp<s<0. 

Thus,  since  Tm  >  g  (see  (4.46)),  one  finds 
(4.70)  <p'(s)  >0  per  s  e  [Tp,T<r] 

where  T p  —  max{TpTm}. 

By  virtue  of  proposition  4.1,  Ts( 0)  belongs  to  the  interval  [Tp,Tj,  if  the  initial  data  are 

consistent  (that  is  if  (4.36)  holds).  Thus,  we  have  v'(Ts(0))  >  0,  and  the  equality  holds  if 
and  only  if  Ts(0)  =  Ta. 

Moreover,  since  Vl(a)  >  0,  y2(s)  >  0  for  Tp  <  s  <  Ta,  we  see  that  the  function  inside  the 
integral  in  (4.68)  is  strictly  positive,  then  Ts(t)  is  increasing,  as  long  as  Ts(t )  <  T(r. 

From  (4.67)  we  deduce  that  Ts(t)  can  not  have  an  asymptotic  value  e<  T<t,  because 
<P2(s)  vanishes  only  for  s  =  T(r.  Consequently,  the  solution  Ts(t )  of  (4.67)  tends 
monotomcally  to  T,,  in  a  finite  or  infinite  time  accordind  to  the  fact  that  the 
integral  m  (4.68)  tends  to  a  finite  or  infinite  value,  respectively  (notice  that  for  s  =  to 
the  denominator  of  the  function  in  the  integral  vanishes,  while  the  numerator  achieves  I 
finite  negative  value).  If  <  oo,  we  expect  a  melting  process  for  t  >  tvy 
If  <oo  =  oo,  that  is 

—  oo 

(4'7I)  » . 

TS(«> 


then  the  solution  of  (S«rap)  +  (C)  +  (A)  (equations  (4.1)-(4.10))  tends  to  stationary  values 
ehich  we  are  going  to  describe. 

The  final  height  of  the  soil  Hx  is  achieved  from  (4.19): 

(4-72)  HOQ  =  H  +  zf-b 


where  zf  =  =  b  + 


kx  g  —  T 

LPi^rf)  =  Fl{T^ 


The  thickness  of  the  frozen  fringe  changes  with  respect  to  time  according  to  the  formula 
(4.1): 
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Ts(t) 

J  kff(v)dV 


(4.73)  b-zF(t)  = 


0 


b. 


LPiVsiT s(P) 
derivating  (4.73),  one  finds: 

(4.74)  iF(t)  =  j  >  0 

(LPl<p3(Ts(t)))2 

Therefore,  the  thickness  of  the  frozen  fringe  decrases,  but  the  isotherm  line  2  =  zF(t) 
does  not  touch  the  base  of  the  lens,  since  6  >  zF(t). 

The  final  thickness  of  the  frozen  fringe  is  given  by 

r, 

J  kj/iiyti 

(4.75)  -y-MU-l'-fespr,- 

Finally,  the  water  flux  qw  satisfies 


feUO  =  o. 

We  conclude  by  remarking  that  the  temperature  Tp  defined  by  (4.67)  is  constant, 
because  of  (4.66);  so,  Tp  =  Tp( 0).  Hence,  condition  (4.22)  is  surely  verified  at  ant  time 
t  e  [0,^),  since  T5(0)  >  Tp  and  Ts(t)  is  increasing. 

We  examine  now  the  general  case  in  which  h  and  g  depend  on  time. 

Let  us  write  explicitily  the  differential  equation  (4.20),  taking  into  account  that  the 
formula  (4.69)  can  be  used  to  compute  the  partial  derivative  : 

01  c 


(4.76) 


where 


-  =  wx{Ts{t),t)Ts{t)  +  w2(Ts(t),t)g(t)  +  w3(Ts(t),t)h(t) 


w  fs  t)= _ kf(H ~ 6)  ( (  kjAfl  1  ^(M)  1  \  ,.s_  ,  , 

1  ’  ip,^3(s7<)(V,(s,<)  +  H-  &)\\&-fy>,-P3(s,P +  ds  ip(s,t)  +  (H -b)J^  '  s>  J 


W2(s,t)  = 


kj£-i) _ 1 

LPi  <P3(sit)Ws-t)  +  H  -b) 


_ AtO-s _ Sil>(s,t)  +  H -b  \ 

(^>3(s,p(t ~p{s~t)  +  77  —  6))  '  ^Lpt  h.QpwLtpi{s,t)) 
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The  function  rf>(s,t)  is  defined  as  in  (4.54): 

(4. 77)  tp(s,t)  =  (H-b)  —  - 

¥>i(W3M 

Obviously,  Wl(s,t)  coincide  with  „'(,)  in  the  examinated  case  h  and  g  constant. 

The  initial  condition  Ts(0)  for  the  differential  equation  (4.76)  is  achieved  from  (4.21), 
which  has  exactly  one  solution  Ts{ 0)  >  Tp  (see  (4.37)),  by  virtue  of  proposition  4  1 

Let  us  assume  the  validity  of  the  following  hypothesis  which  is  absolutely  not  restrictive 
for  the  generality  of  our  analysis: 

(4.78)  sup  h(t)  =  h<  oo 
t  >  o 

We  define 


(4.79)  Tp  =  Tp(h) 

where  Tp  satisfies  the  equation  (4.67). 


PROPOSITION  4.2.  Assume  that  (4.78)  holds.  If  the  functions  k„(s),  Kl(s), 
X2(s)eC  (-CO.0],  4(1),  »(l)«C*(0,oo)  and  if  the  initial  data  fulfil  condition  (4  36)  then 
there  exists  an  interval  of  time  [.,<,)  uthere  (4.4.76)  has  a  unique  solution  T M)  in 
iT p{h),T a],  where  Tp(h)  is  defined  by  (4.79)  and  T  by  (4  34) 

Dim . 

From  (4.36)  we  have  (see  remark  4.2): 


^>T5(0)>Tp(A(0))>rp(^). 

We  introduce  now  the  function 

(4-80)  =  -Tr^J  (£$r +  +  r3(«.04(o). 

Computing  «(F(M)/».,  one  can  easily  check  that  €  C'[rp(S),TJ,  on  the  ground  of 

the  assumed  hypotheses  for  kff(s),  Kl(. s),  K2(s),  h(t)  e  g(t)  and  of  the  estimates 
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(4.81) 


Pib 


f  zu_ 

J  «i 


kff(v) 
0 7 ) 


TJh) 


dr) 


kff(v) 


PwK o  TP(A) 

M  -  J  kff(v)dr] 


p-b 

<  <Pi («,<)<  r-g- 

r  0 


/  *'(„) 

o _ 

T* 

~  J  kff(v)drl 


■dr) 


(4.82) 


Tp(h)  ) 

kuh(t)-  j  kff(rl)drl 
0  / 


By  virtue  of  the  local  existence  and  uniqueness  theorem  for  ordinary  differential 
equations,  we  conclude  that  there  exists  a  finite  time  tj  >  0  such  that  (4.76)  has  exactly 
one  solution  Ts(t)  with  the  initial  datum  Ts( 0)  and  whose  graphic  is  contained  in  the 
compact  set  [0,fy]  x  [Tp(/i),T  J.  □ 

Once  (4.76)  has  been  integrated,  we  can  compute  the  boundaries  zF{t)  and  zT(t)  from 
(4.12)  and  (4.19),  respectively.  The  water  flux  is  given  by 


<7  JO  = 


^2(T5(0) 

PwViiTsiW 


The  height  of  the  soil  T(t)  verifies  also  the  formula 


(4.83) 


zT(t)  —  H  ~ 


f  <P2(Ts(t)) 
J  <Pi(Ts(t),t) 
0 


dr. 


Finally,  the  water  pressure  pw  is  given  by  (1.28a). 

We  are  going  now  to  check  the  sign  of  the  coefficients  i  =  1,2,3,  in  equation  (4.76). 


PROPOSITION 4.3.  We  have: 


(4.84)  W^sJ)  >  0  for  Tp(h(t))  <  s  <  0,  t  >  0 

(4.85)  W2(s,t)  <0  for  5  <  Ta,  t  >  0 

(4.86)  W3(s,t)  <  0  for  g(t)  <s<Ta 

where  Tp{h(t))  fulfils  (4.25)  at  each  time  t. 

Dim. 

(4.85)  comes  from  (4.18),  (4.40)  and  (4.41). 
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Moreover,  (4.84)  is  proved  by  remarking  that 


dip(s,t) 

ds 


=  «'(*) 


wheref(s)  is  defined  by  (4.54)  and  by  recalling  (4.18),  (4.55),  (4.57). 
Let  us  show  now  (4.86). 

If  we  derivate  (4.15)  with  respect  to  s ,  we  find: 


(4.87) 


fyiM  Pibkff(s ) 

d*  PWK o 


Since  Kx{s)  is  increasing  and  K^O)  =  K0,  it  follows  that  <  0  for  a  <  0,  t  >  0. 

Therefore: 


¥>i(s,t)  <  ^(0,0  =  for  s  <  0,  t  >  0 
Pwn  0 

and  consequently 

bLPi  KoPwLVl{s,t)  -  ^t\ bP,  -  J  >  0  for  s  <  0,  <  >  0 

and  (4.86)  is  thus  proved  □ 

PROPOSITION  4.4.  Assume  that  Ts(t)  is  solution  of  (4.76),  Ts{t)  e  [Tp(h),Tv];  then,  for 
each  time  t  e  [0,^]  we  have  Ts(t)  >  g(t). 

Dim .  Call 


*.  =  *«/{<€  [<>,<,):  T5(<)  <*(i)}. 

By  (4‘36)  we  see  that  TsQ)  >m  We  ta  >  0.  By  the  definition  of  inf  we  have 
Ts(tJ  =  9(ta)  and  Ts(t)  >  g(t )  whenever  t  €  [0,<o).  it  follows  that 
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(4.88)  Ts(ta)  <  g(ta). 


By  the  definition  of  the  functions  W2  and  W3  one  easily  finds: 


(Ts(ta),ta),  WaW  =  0. 


Therefore: 

(4's9)  “  =  -  ~  «'•>) 

Suppose,  contrary  to  our  claim,  that  Ts(ta)  <  Ta.  In  that  case,  we  would  have  that  the 
left-hand  quantity  in  (4.89)  is  strictly  positive  (see  (4.40)  and  (4.41)),  while  the  right- 
hand  one  is  non  positive,  due  to  (4.88)  and  to  proposition  4.3.  Thus,  we  would  obtain  a 
contradiction.  If  Ts(ta)  =  Tff,  the  left-hand  side  of  (4.89)  vanishes,  but  we  can  conclude 
as  in  the  previous  case  by  considering  times  t  close  to  *a?  t  <  ta  (in  (4.88)  holds  the  strict 
inequality,  for  such  as  times).  □ 

Let  us  consider  now  the  maximal  interval  where  the  solution  Ts(t)  of  (4.76)  is  defined, 
in  the  sense  of  proposition  4.2;  we  keep  on  calling  tf  the  right  boundary  of  the  time 
interval. 

PROPOSITION  4.5.  Under  the  same  assumptions  as  in  proposition  4.2,  the  solution 
Ts(t)  o/(4.76)  satisfies  one  and  only  one  of  the  following  possibilities : 

1)  Ts(tf)  =  T,,  and  Ts(t)  >  Tp(h(t))  Vf  6  [0,^]; 

2)  3  t1e  (0 such  that  <  Tp(h(tf))\ 

3)  Tp(t)  <  T s(t)  <  T g  V<  >  0. 

Dim.  It  suffices  to  take  into  account  propositions  4.3,  4.4  and  recalling  that  the 
maximal  solution  can  not  be  contained  in  any  compact  set  which  is  properly  contained 
in  [O,^] x[Tp(h),Ta}.  □ 

Let  us  comment  the  three  possibilities. 
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In  the  case  1),  the  solution  Ts(t)  reaches  the  value  Ta  and  the  water  flux  qw  vanishes  at 
that  temperature.  For  t  >  tf  condition  (4.24)  is  violated  and  a  melting  process  occurs. 

In  the  case  2),  by  the  continuity  of  Ts(t)  and  by  the  estimate  Tp(h)  <Tp(h(0))  <TS(0), 
there  exists  a  time  t2  such  that  T s(t2)  =  T p(h(t2)).  The  solution  Ts(t)  cal  be  accepted 
only  for  t  <  t2,  since  for  t  >  t2  the  water  pressure  gradient  at  the  base  of  the  ice  lens 

becomes  negative.  It  is  to  be  expected  that  for  t  >  t2  a  process  frost  penetration  will  take 
place. 

The  case  3)  occurs  when  tf  =  oo  and  the  solution  Ts(t)  remains  bounded  by  denoted 
values.  We  have  already  met  with  this  possibility,  when  the  temperatures  h  and  g  are 
constant  and  the  integral  defined  in  (4.68)  verifies  (4.71).  The  solution  Ts(t)  tends 
monotonically  to  T  a.  However,  other  asymptotic  values  for  the  solution  Ts(t)  are 
possible,  m  the  case  tf  =  oo:  as  an  istance,  we  can  choose  in  (4.76)  h  constant  (in  that 
case  W1  and  W^2  don’t  depend  on  time  explicitily)  and  g( 0)  such  that  (4.35)  holds.  Once 
the  initial  value  Ts(0)  is  calculated  by  (4.21),  we  define 


(4.90) 


9(t)  =  g(  0)  - 


^i(2V0)WVT5(0)r 


It  is  easily  seen  that  the  solution  of  (4.76)  is  Ts(t)  =  T5(0). 

Let  us  now  deal  with  the  following  question:  how  the  boundary  temperatures  h(t),  g(t) 
has  to  be  chosen  in  order  to  discriminate  the  three  possibilities  just  described? 

For  this  purpose,  we  state  the  following  results  (propositions  4.6  and  4.7),  which  allow 
us  to  answer  the  question  just  introduced  at  least  in  special  cases. 


PROPOSITION  4.6.  Let  assumptions  stated  in  proposition  4.2  hold.  Then,  if  h(t)  and 
g{t)  are  non  decreasing,  then  the  solution  Ts{t)  of  (4.4.76)  is  strictly  increasing  for 
0  <  t  <  tf  and  it  reaches  in  a  finite  time  or  asymptotically  the  temperature  T0. 

If  h(t)  and  g(t)  are  non  increasing  and 


+  (h(t))2  >  o, 


0  <  t  <  tf, 


then  Ts(t)  <  Ta. 


Dim.  Assume  that 
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(4.91)  g(t)  >  0,  h(i)  >  0. 

By  virtue  of  (4.91),  the  temperature  Tp(h(t))  is  non  increasing  (see  remark  4.2). 

On  the  other  hand,  by  propositions  4.3  and  4.4  and  by  (4.40),  (4.41),  we  see  that  the 
function  W  defined  by  (4.80)  verifies  at  each  time  t  such  that  the  solution  Ts(t)  of  (4.76) 
exists 

(4.92)  W(Ts(t),t)  >0,  0  <  t  <  tf 

it  follows  that  Ts(t)  is  non  decreasing.  Since  T5(0)  >  Tp(h( 0))  and  Tp(h(t))  is  non 
increasing,  we  obtain 

(4.93)  Ts(t)  >  Tp(h(t)). 

Thus,  it  is  not  possible  that  the  second  possibility  of  proposition  4.5  occurs. 
Furthermore,  from  equation  (4.76)  we  deduce  that  the  solution  Ts(t)  can  not  have 
asymptotic  values  lower  than  the  temperature  Ta:  it  may  be  concluded  that  the 
temperature  Ts(t)  reaches  (in  a  finite  or  infinite  time)  the  value  Ta. 

We  show  now  that  Ts(t)  is  strictly  increasing  for  0  <t<tp 

The  equality  in  (4.92)  holds  if  and  only  if  Ts(r)  =  T g{r)  =  0,  h(r)  =  0,  for  some  r  e  [0,^]. 
Since  Ts(t)  is  the  maximal  solution,  it  must  be  r  —  tp  hence 

Ts(t)  =  W(Ts(t),t)  >0,  for  0  <  t  <  ip 

Let  us  prove  now  that,  keeping  the  assumption  (4.91),  the  solution  Ts(t)  has  the 
following  property: 

Ts(t)>T°s(t),  0  <t<t} 

where  T%(t)  is  the  maximal  solution  of  (4.76),  defined  in  [0,ty),  obtained  by  taking  the 
boundary  temperatures  as 

(4.94)  h°(t)  =  h(0),g°(t)=g(0). 
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The  solutions  Ts(t )  and  T$(t)  satisfy,  respectively: 


Ts(t)  =  W(Ts(t),t) 
T°s(t)  =  W0(T°s(t)) 


where  W(s,t )  is  defined  by  (4.80)  while  W0(s)  =  -y2(5)/^(s)^(s)  (see  (4.67)). 

For  the  sake  of  clearness,  we  recall  that  the  functions  <p,  y,  with  just  one  argument  are 
related  with  the  examinated  case  h  and  g  constant. 

We  remark  that  the  initial  datum,  that  is  the  solution  of  (4.21),  is  the  same  for  both  the 
problems  (so  Ts( 0)  =  by  virtue  of  (4.94). 

Let  us  integrate  both  the  differential  equations  in  a  interval  where  both  the  solution 
exist: 


Ts(t)  T°s(t ) 

J  w&ify  =  J  = t- 

TS(  0)  Ts(  0) 


and  write  explicitly  W  and  W0: 


w{y’t)=  ~w^i{  +^m*)+w^tm) 


w°(y)  =  - 


i  ^(y) 
wm  <pM 


By  assumption  (4.91),  proposition  4.3  and  (4.93),  we  get 


(4.95)  -  W2(y,t)g(t)  -  W3(y,t)h(t)  >  0 


On  the  other  hand,  from  the  definition  of  Wit  *  =  1,2,3  and  keeping  in  mind  (4.91)  once 

again,  it  follows  immediately  that,  for  each  time  such  that  both  the  solutions  of  the  two 
problems  exist,  it  is 

(4.96)  2^<„,^f!<p2<0 
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We  conclude  that  W(y,t)  >  W°(y)\  therefore  Ts{t)  >  T%(t). 

Roughly  speaking,  we  may  say  that,  under  the  same  initial  condition  Ts(0),  the 
increasing  boundary  temperatures  ”  bring”  the  solution  Ts(t)  to  the  value  T a  in  a  shorter 
time  with  respect  to  the  case  when  the  boundary  temperatures  are  settled  at  the  values 
A(0),  <?(0). 

PROPOSITION  4.7.  If  h(t)  and  g(i)  satisfy  the  conditions 
(4.96)  g(t)  <  0,  h(t)  <  0,  ( g(t ))2  +  (h(t))2  >0,  0  <  t, 

then  Ts(t )  <  T c. 

Dim.  Assume  that  (4.96)  holds.  If  there  existed  a  re[0,*y]  with  the  property  Ts(r)-T^ 
we  would  have  T5(r)  >  0.  By  propositions  4.3,  4.4  and  assumption  (4.96),  we  would  have 
W(Ts(t),t)  <  0,  where  W  is  defined  by  (4.80)  and  we  would  obtain  a  contradiction.  □ 
the  example  (4.90)  shows  that,  even  if  the  temperatures  h  and  g  are  non  increasing  (in 
the  example  the  temperature  g  is  really  decreasing  in  a  linear  way),  the  temperature 
Ts(t)  is  not  necessarily  decreasing:  the  term  in  the  left-hand  side  if  equality  (4.76)  is 
positive,  while  in  (4.95)  it  has  the  opposite  sign.  In  qualitative  terms,  we  can  say  that,  if 
g(t)  decreases  quickly,  the  temperature  Ts(t)  reaches  (by  decreasing)  the  value  Tp(h(r)) 
for  some  re(0,iy)  only  if  g(t)  decreases  rapidly  enough.  Moreover,  when  h(t)  decreases, 
the  temperature  Tp(h(t))  increases. 

Let  us  give  now  the  following  example. 

We  choose  h(t)  =  /i0,  where  h0  is  a  positive  constant  and  ^(0)  is  taken  such  that  (4.36) 
holds.  By  means  of  (4.21)  we  calculate  T5(0)  and  we  define: 

f(t)  =  at  +  Ts( 0),  a  <  0. 

Since  h  is  constant,  the  temperature  Tp  which  satisfies  (4.25)  is  also  constant; 
moreover,  the  functions  W2  and  ^  don’t  depend  explicitly  on  time. 

It  is  a  simple  matter  to  find  a  and  g(t)  such  that 

(4.97)  /(,) =«>»-(/«),  o = -  + w) 
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Inequality  (4.97)  means  essentially  that  g(t )  must  decrease  rapidly  enough. 

The  solution  of  (4.76)  that  is  achieved  by  imposing  the  boundary  temperatures  h0  and 
s(<)  (where  g(t)  verifies  (4.97))  keeps  below  the  straight  line  /(<).  The  temperature  TM) 

g  necessarily  down  the  vallue  Tp(h0)  and  the  water  pressure  gradient  at  the  base  of 
the  lens  becomes  negative. 

Owing  to  condition  (4.22),  the  solution  describing  the  process  of  lens  formation  can  be 
accepted  only  up  to  the  time  t*  when  Ts(t*)  =  Tp(h0). 

4.2  Frost  penetration 


In  this  section  we  wiil  look  for  solutions  of  (Stmp)  +  (Q  +  (A)  (eqq.  (4.1)-(4.1))  which 
describe  the  process  of  frost  penetration ;  since  it  must  holds  (see  par.  4.1.7) 

(4.98)  is(t)  <  o, 


from  (4.4)  and  (4.7)  we  get 

<4'99)  zvfflhm 

From  (4.1),  (4.4)  and  (4.99)  we  deduce  the  equation  for  the  temperature  Ts(t): 


T s(t)  k  (n\  ^2 (T s(t)) ,  ,  v 

(4.100)  a  +  k_MMLsW)  ,  [  M 

v  y  ~  .  7t  f+\\  +  /  - - - ,  - - -An  — 


Ko  kfj{Ts{t)) 


KM 


-dr)  =  0. 


Equation  (4.100)  is  equivalent  to 


(4.101)  V>2(r5(t))  pwL  kff(TSg(t))  s{t),t), 

(where  <puV>2  and  <p3  are  defined  by  (4.15)-(4.17)),  or  even  to 
(4.1.02)  r(Ts(t))  =  kuh(t ) 
where  r(s)  is  defined  by  (4.27). 

The  solution  of  (4.102)  is  given  by  the  temperature  Tp(h(t ))  defined  by  (4.25).  Therefore, 
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we  assume  that  the  temperature  Ta  verifying  (4.34)  exists,  otherwise  (4.102)  has  no 
solution  (see  lemma  (4.25)). 

Formally,  once  (4.102)  has  been  solved,  the  freezing  temperature  Ts,  which  depends  on 
time  t  only  through  the  boundary  temperature  h(t).  is  a  known  function. 

Writing  qw(t)  and  zF(t)  in  terms  of  zs(t)  and  zT(t)  by  means  of  the  formulas  (4.1),  (4.4) 
and  (4.99),  one  finds: 

(4.103)  qjt)  =  -  ^-A-^^Lp{ip3(Ts(t),t)  j^-t y 

(4.104)  zF(t)  =  y(t)zs(t) 


(4.105)  y(t) 


kuh(t) 

Ts(t) 

kuh(t)~  J  kff(tl)dri 


we  have  0  <  7  <  1,  since  if  Ts(t )  is  the  solution  (4.100),  it  must  be  Ts(t)  <  0. 

Substituting  (4.103)  and  (4.104)  in  (4.2)  and  in  (4.4.3),  we  get  the  following  set  of 
ordinary  differential  equations,  where  the  unknown  quantities  are  the  boundaries  zs(t) 
and  zT(t): 

1  fL9(t)~Ts(t)  PwLK2(Ts(t))  +  kff(Ts(t))ir  . .  j] 

s{t)  ~  (1  -  us)epwL  |*/zr(<)  -  zs(t)  kff(Ts(t))  *LPi<p3{Ts{t),t) 

;  a  (t)-kM~PJ  9(t)-Ts(t )  bLPi<p3(Ts(t),t)  f 

ZT(t)  ~  —pj  |  PiPwLK2(Ts{t))  + 


+  (Pi-pJkff(Ts(t)) 


The  initial  conditions  are  given  by  (4.6): 


rs(0)  =  6,  zT( 0)  =  H  >  b. 


By  means  of  the  definiton  of  the  variables 
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C(t)  =  _  {p'PwLK^TS^))  +  (Pj  -  Pw)kff(TS(t))}(  1  -  vs)e  -  Pi{kff(Ts(t))  +  PwLK2{Ts{t))] 

~  PiPw^~vs)kff(Ts(t))  ~  'x 

x  bPi<P3(T s(<),<)  =  1  -  us)e  - 1)  +  bV3{T s(t),t) 

D(t)  =  p)  _  kf(g(t)  -  Ts(t))(Pi  -PJ  _ 

PiPwLe{l-vs)  PtPwL{^  -  t's) 

The  initial  condition  for  SF  are 

(4.106)  x(0)  =  b,y(0)  =  H-b 

Let  us  now  discuss  the  signs  of  the  just  defined  coefficients.  From  (4.40)  we  obtain: 

(4.107)  A(t)  >  0,  C{t)  <0  for  s  <  0,  <  >  0. 

Moreover,  noticing  that  ((1  -us)c{Pi-PJ  -Pt)<  0  (this  is  true  even  in  the  case  Pw<Pi), 
we  have: 
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(4.108)  B(t)  <  0,  D(t)  >  0  for  g(t)  <  Ts(t). 


Analogously  eith  the  case  of  lens  formation,  we  will  deal  first  the  problem  to  check  in 
which  cases  the  initial  data  are  consistent  with  the  prescribed  conditions  (4.10),  (4.11), 
(4.98). 

Under  the  assumption  (4.34),  which  is  assumed  to  be  satisfied,  we  know  that  the 
equation  (4.102)  has  one  solution  for  each  fixed  time  t . 

Let  us  call  to  shorten  notation  Tp(h( 0))  =  Tp(0)  and  state  the  following 

PROPOSITION  4.8.  The  initial  temperature  T5(0)  =  Tp(Q),  solution  of  (4.102)  calculated 
for  t  —  0  is  consistent  with  (4.10),  (4.11): 


0  <  zF(t)  <  zs(t)  <  zT(t ),  qw(t)  >  0 


and  (4.4.98)  if  and  only  if 


(4.109)  kfg( 0)  <  k,Tp(0)  +  (H-b)LP{<p3(Tp(0))  ^1  + 

Dim. 

Assume  that  (4.109)  holds. 

From  (SF)  we  get 


(4.110) 


.m_  i  fkjjm-Tsm 

[  )~  &  +n-b  -(i-KTsm>pw{  L(S-b) 


PwLK2(Tsm+k„(Tsm 

epwkff(Tsm 


Pi^Tsi0)’0) 


<0. 


Hence,  condition  (4.98)  is  fulfilled  for  t  =  0. 

Condition  (4.10)  is  obviously  satisfied  for  t  =  0  (see  (4.104)). 

Furthermore,  the  solution  of  (4.102)  must  benegative;  thus,  from  (4.103)  and  from 
(4.40)  we  see  that  ^(0)  >  0,  that  is  (4.11)  for  t  =  0. 

On  the  other  hand,  if  (4.109)  were  not  true,  from  (4.110)  we  would  find  %  (0)  >  0, 
contrary  to  condition  (4.98).  □ 


Keeping  in  mind  (4.98)  and  (4.10),  we  see  that  the  solution  (x(t),  y(t))  of  (SF)  must 
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verify: 


(4.111)  x(t)  <  0, 

(4.112)  x(t)  >  0,  y{t)  >  0. 


Condition  (4.11)  is  certainly  fulfilled  if  (4.11)  holds  (see  (4.40)  and  (4.103)). 

Generalizing  (4.110)  to  any  time  t  and  taking  into  account  the  equation  (4.102),  we  see 
that  (4.111)  is  equivalent  to 

(4.113)  k„[t)  <  wm  +  ( 1  +  I 

V  k  f  M  pW)) ) 

where  Tp(h{t))  =  Ts(t), 


or  to 


(4.114) 


y(t)  B(t) 

x(t)  A(ty 


If  Tg(t),  solution  of  (4.102),  satisfies  (4.113),  it  must  be 
(4.115)  Ts(t)>ff(t). 

Let  us  now  calculate  the  quantity  AD-BC;  an  easy  computation  shows  that: 

(4.H6)  a(wo  -  mm = iPj 

^  Ts{t)  is  solution  of  (4.102),  we  get  from  (4.115) 

(4.117)  A(t)D(t)-B(t)C(t)>0. 


In  order  to  discuss  the  problem  (Sf),  let  us  investigate  first  the  easier  case  a  and  , 
constant. 
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PROPOSITION  4.9.  Assume  h>  0  and  g<  0  costant  and  let  be  (4.109)  verified .  Then , 
there  exists  a  time  t0<  oo  such  that  (S^)  has  exactly  one  solution  (x(t),y(t))  in  the  interval 
*e[0;*o),  50  that  (4.111),  (4.112)  are  satisfied  for  *e[0;*o);  while  x(t0)  =  0. 

Dim . 

The  system  (SF),  in  the  present  case  of  constant  boundary  temperatures,  is  autonomous 
(A^BjC^D  are  constant),  since  the  temperature  Ts  solving  (4.102)  is  constant. 

Owing  to  condition  (4.112),  we  are  interested  just  in  positive  solutions  fn(/),i/(*)). 

By  the  theorem  of  existence  and  uniqueness  for  ordinary  differential  equations,  we  can 
say  that  for  each  pair  of  initial  values  z(0)  >  0,  y(0)  >  0  there  exists  a  unique  solution  of 
(SF)  at  least  locally. 

Consider  now  the  projections  of  the  solutions  fr(2), ?/(£))  on  the  quarter  of  plane 
Q  =  {x  >  0,y  >  0}. 

Let  us  write  the  orbits  in  Q  in  the  form  y  =  y(x ).  We  find  the  differential  equation 

where  the  apex  denotes  the  derivative  with  respect  to  x. 

Conditions  (4.111),  (4.112)  impose  that  the  starting  point  (x(0),y(0))  has  to  lye  in  Q  and 
that  (4.114)  evaluated  for  t  —  0: 

(4.119)  ^<-f 
must  hold. 

Integrating  (4.117)  and  taking  into  account  of  the  initial  conditions  (4.106),  one  finds 
the  following  formula: 

(4.120)  -  ±ln j  -  Au\t )  +  (C -  B)u(t)  +  D | -  cx In  f  ^  =/n^|+c4 


where  we  set 


c2  —  B  —  C  —  $AD  +  (5  -  C)2, 


c3  =  5-C  +  a[4AD  +  CB^C)2, 

C4  =  -  -  MO)  +  (c  -  5)11(0)  +  D)~  Clln^;cl\ 

with  u(0)  =  -  ~  h. 

We  remark  that  c3  >  0  and  that  (4.117)  and  (4.119)  yeld 
c2  25,  —  Au^( 0)  +  ( C  —  5)i4(0)  +  D  >  0. 

Therefore,  the  constant  c4  is  well  defined. 

We  see  that  y'  =  0  if  and  only  if  the  point  (x,y)  belongs  the  half-straight  line 
(4.121)  s0  =  {(x,y)eQ  :  y= -Dx/C}; 


moreover,  y‘  =  00  if  and  only  if  (x,y)  lies  on  the  half-straight  line 
(4.122)  &<„  =  {(x,y)eQ  :  y  = -Bx/A}. 

By  virtue  of  (4.117),  the  slope  of  s0  is  greater  that  the  slope  of  s^. 

If  we  look  for  half-straight  lines  y  =  mx  in  Q  which  are  orbits  for  the  system  (SF),  we 
easily  see  that  m  must  verify  Am2  +  (5  -  C)m  -  D .=  0.  Such  as  equation  has  two  solutions: 

mi=~l i<°,  m2=~ |>0. 

Moreover:  -j<m2< 

Thus,  inside  the  angle  bounded  by  s0  and  by  Soo  there  exists  one  straight  line-orbit. 
Define  the  three  angles 

ri  =  {(*>2/)  e  Q  :  0  <  y  <  -  Bx/A} 
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r2  =  {(^J/)€Q  :  —  Bxj  A<y  <  -Dx/C} 


r3  =  {(x>2/)  6  Q  :  y  >  -Dx/C}. 


It  is  easily  seen  that  yf  >  0  if  and  only  if  (x,y)  e  T2.  Moreover,  x  <  0  if  and  only  if  (x,y)  e  rx, 
y  >  0  if  and  only  if  (x,y)  e  Tt  U  r2. 

Owing  to  condition  (4.111),  we  have  to  consider  only  the  orbits  in  the  angle  Tv  where 
x  <  0,  y  >  0  hold.  Choose  an  initial  point  P0  =  (z(0),t/(0))  e  and  follow  the  orbit  starting 
from  P0.  As  long  as  the  orbit  remains  in  r2,  we  have  x(t)  <  0,  y(t)  >  0.  The  orbit  meets 
necessarily  the  half-straight  line  y  =  -  Bxj A  in  a  finite  time  t  =  <0;  at  that  time  *0,  the 
boundary  zs  is  at  the  height: 


^sC^o)  —  be 


AD-BC 


(2B  —  c2\ 
\c3-2B) 


>0. 


We  have  £(*0)  =  0  and  for  t  >  t0  it  is  x(t)  >  0  and  condition  (4.111)  is  violated.  We  will 
check  later  if  for  t  =  t0  a  process  of  lens  formation  takes  place.  □ 

Once  the  system  ( SF )  has  been  integrated,  the  water  flux  qw  can  be  computed  by 
(4.103): 

The  function  7  defined  by  (4.105)  is  constant: 


(4.123)  zF(t)  =  yx(t) 


with 


0  <  7  <  1. 


From  (4.117)  we  deduce  y(t)  >  0  for  <c[0,to);  taking  also  into  account  (4.123),  we  find 
that  the  condition  (4.10)  is  also  verified  in  the  interval  t  e  [0,*o). 

We  remark  that  the  thickness  of  the  frozen  soil  y(t)  —  zT{t)  -  zs(t)  is  increasing  for 
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0  <  t  <  t  by  virtue  of  (3.48). 

The  height  of  the  soil  at  the  time  t  =  is  given  by 
zT(f  o)  =  v(t  o)  +  *(<0)  =  (i  -  jjzsCto)’ 

Let  us  consider  now  the  system  (SF)  in  the  general  case  when  the  boundary 
temperatures  h  and  g  depend  on  time.  As  in  the  previous  case,  we  get  the  temperature 
Ts  from  (4.102),  but  m  the  present  case  Ts  depend  on  time:  in  particular,  if  h(t) 
increases  (decreases),  the  temperature  Ts  decreases  (increases). 

Let  us  exammate  once  again  the  quarter  of  plane  Q  and  especially  the  angle  r,f<) 
defined  by  1  ’ 


ri(<)  =  {(«,y)e<3  ••  0  <y<  -  B(t)x/A(t )}; 


Likewise  the  case  h  and  g  constant,  (4.111)  and  (4.112)  evaluated  for  t  =  0  impose 
po  =  (^(0),y(0))  eT1(0),  that  is 


(4.124) 

V  '  b  ^(0)' 


Assuming  that  the  prescribed  functions  kff(s),  K^s),  k2(s),  h(t)  and  g(t)  have  the  same 
regularity  as  stated  in  proposition  4.3  and  applying  the  theorem  of  existence  and 
umqueness  for  ordinary  differential  equations,  we  have  that  one  solution  (s(<),y(<))  of 
(SF)  starting  from  P0  is  defined.  Consider  the  projection  of  the  solution  on  Q  (obviously, 
m  the  present  case  the  orbits  are  no  longer  of  the  form  (4.119)).  As  long  as  the 
projection  remains  in  the  angle  r^t),  whose  width  depends  on  time,  the  conditions 

(4.10) ,  (4.11)  and  (4.98)  are  satisfied.  Indeed,  is(t)  <  0  if  and  only  if  y/x<  -  B{t)/A{t ); 

(4.11)  is  trivially  verified  (see  (4.103)).  By  (4.117)  and  by  (4.104)  (4.10)  is  also  fulfilled.’ 
Let  us  write  (4.110)  in  the  following  way: 


(4.125) 


y(0  £(<) 

*(<)  A(t)  - 


bm~Ts(t)) 


Lptbtf3(Ts(t).1) 


1  +  Pwp 


F(h(t),g(t)) 


According  to  what  we  said  above,  the  process  of  frost  penetration  stops  when  and  only 
when  in  (4.125)  holds  the  equality  for  some  time  r: 
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(4.126) 


y(r)_  b(t) 
x(t)  A(t) 


if  and  only  if  zs(T)  — 


From  a  geometrical  point  of  view,  (4.126)  holds  when  the  projection  of  the  solution 
(y(t),x(t))  matches  the  straight  line  y—  -B(t)x/A(t),  which  is  the  upper  side  of  the  angle 

r,(*)- 

The  ratio  y(t)/x(t)  is  certainly  increasing  if  ( x(t),y(t ))  is  solution  of  (SF)  with 

^-(xiowo^r^o). 

Nevertheless,  the  behaviour  of  the  function  F  with  respect  to  time  is  related  to  the 
profile  of  the  boundary  temperatures  h(t)  and  g(t). 

We  state  the  following  result. 


PROPOSITION  AID.  Assume  that  (4.124)  is  verified.  If  F(h(t),g(t))  is  finite  for  t>  0, 
then  there  exists  a  r  >  0  such  that  (4.126)  holds. 

Dim. 

By  hypothesis  we  have 


(4.127)  sup  2<oo 

t  >  0  A\t) 

The  projections  on  Q  of  the  solutions  (x(t), y(t))  of  SF  satisfy  the  equation 


(4.128) 


y'O)  = 


Cft)y  +  Dftyc 
A(t)y  +  B(t)x' 


The  isocline  straight  line 


soo(t)  =  {(*>y)  e  <3  :  y  =  -  B{t)x/ 4(<)} 

has  a  positive  finite  slope,  depending  on  time;  by  the  assumption  (4.127),  the  slope  of 
soo(t)  *s  finite. 

Condition  (4.124)  assures  that  (r(0),t/(0))  €1^(0).  Moreover,  from  (4.128)  we  see  that  the 
projection  y  =  y(x)  is  decreasing  with  respect  to  x  as  long  as  (x(t),y(t))  remains  in  r^*). 
On  the  other  hand,  y  —  y{x)  can  not  accumulate  in  any  point  in  r2(<),  since  thera  are  no 
singular  points,  owing  to  (4.16).  Thus,  the  projection  of  the  solution  must  reach  in  a 
finite  time  t  —  r  the  straight  line  ^(r)  and  ir(r)  =  zs(r)  =  0.  □ 
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Remark  4.3.  The  result  obtained  by  proposition  4.10  can  be  applied  in  the  following 
special  cases: 

0  h(t)  and  g(t)  are  non  decreasing  (in  particular  constant):  actually,  the  temperature 
Ts(t)  is  in  that  case  non  increasing  and  the  function  F  is  bounded  by 


F(h(i), g(t))  < 


<?(0)  -  Ts(0) 
P3(°) 


t  >  0; 


n)  g(t )  is  lower  bounded:  indeed,  if  g(t)  >  -gQ  >  _oo,  t  >  0,  we  have 


mt), g(t))  < - Jllh- 


t  >  0; 


u 

ku  /  ktfWl 
0 


dF(h(t),g(t)) 

111  >  di  ^  0:  m  that  case  F(h(t),g(t))  <  F(h(0),g(0)). 


Once  (SF)  has  been  solved,  the  water  flux  gjt)  and  the  isotherm  z  =  zF(t)  can  be 
achieved  by  means  of  (4.103)  and  (4.104),  respectively. 

Arguing  as  in  the  case  h  and  g  constant,  one  can  easily  check  that  the  conditions  (4.10) 
(4.11)  and  (4.98)  are  satisfied  for  0  <  t  <  r. 

One  may  wonder  how  the  boundary  temperatures  h  and  has  to  be  chosen  in  order  to 
have  a  process  of  frost  penetration  for  any  time  t  >  0.  Remarking  that  if  is  bounded 
mt),g(t))  is  bounded,  too,  for  any  value  of  h(t),  we  may  say  that  it  is  necessary  that 

inf  g(t )  =  -  oo. 
t>0 

In  qualitative  terms,  the  just  written  condition  corresponds  to  a  rapid  freezing  process, 
as  we  expect  in  a  process  of  solely  frost  penetration. 

With  regard  to  that  possibility,  let  us  give  the  following  example.  Consider,  for  the  sake 
o  convenience,  t  =  1  as  initial  time  and  precribe  the  boundary  temperatures 


h(t)  =  h0>  0 


where 


s(t)  =  Tp(h0)- 


UfLC-AD(2  n,Db[l 
tf,  2bB  "  l>+~i 
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tl=(l -■?)*<■-*  >0. 


Since  the  temperature  h  is  constant,  we  have  that  T  ,  A  and  C  are  constant  (see  (4.101) 
and  the  definition  of  the  coefficients  A,  C).  Furthermore,  the  ratio  D/B  does  not  depend 
on  g(t),  hence  neither  ( BC-AD)/B  does.  We  find  g(t)  <  0  by  virtue  of  (4.117).  Condition 
(4.119)  is  verified  anyhow  the  constant  h0  >  0  is  chosen. 

The  system  SF  has  the  form: 


(SF) 


ti_A  ,WKI 

\x  —  x  +  y 

c  k2(g(t)  -  T  ) 

y=%+ - v — — 


*(!)  =  &,  y(l)  =  H~b, 


where 


kf{{l-vs)e{Pi-pw)-pt) 

k2  =  - - m - \ - <  0 

PiPwLeO-  ~  J's) 


It  easily  seen  that  the  solution  of  SF  is 

m  =  bp  y(t)  =  +  ^{j- 0+  H - <>1- 


The  boundary  x{t)  =  z5(<)  tends  asymptotically  to  the  base  of  the  soil  z  =  0;  the  velocity 
x  never  vanishes  for  all  t  >  1:  therefore,  we  get  a  process  of  frost  penetration  for  each 
time  /  >  1  and  the  occurrence  of  the  formation  of  an  ice  lens  is  precluded.  The  boundary 
y(t)  is  in  that  example  positive  and  increasing,  since  from  (4.117)  we  get 


•  / 1\  k-tC  ~~  knA  ktyb 

»W  =  JTE7i-!-A>0- 
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4-3  Transition  from  one  process  to  the  other 


Putting  together  the  results  we  stated  in  propositions  4.1  and  4.8,  we  are  able  to  foresee, 
on  the  ground  of  the  knowledge  of  the  initial  boundary  temperatures  h( 0)  e  <?(0)  and  of 
the  properties  of  the  soil  KfT),  K2(T)  and  kff(T),  whether  a  process  of  lens  formation, 
or  frost  penetration  or  melting  will  occur.  As  a  matter  of  fact,  assuming  that  the 
temperature  Ta  defined  by  (4.34)  and  solving  the  equation  (4.26)  evaluated  for  t  =  0  we 
find  the  temperature  Tp(h( 0)).  At  this  point,  if  the  condition 


kf9(0)  >  kfT 0  +  {H  -  b)Lpitp3(T(r,0), 


holds,  where  H  is  the  known  initial  heigiht  of  the  soil,  6  the  known  initial  position  of  the 

freezing  front  zs  and  y>3  is  defined  by  (4.17),  then  the  temperatures  h  and  g  do  not  give 
rise  to  any  freezing  process. 

On  the  contrary,  if  Tp  =  Tp(h( 0))  is  such  that 


(4.129)  kfTp  +  (H  b)LPi<pz(T p,0)(l  +  Pw£j^J^j  <  <kfTa  +  (E  -  b)LPi<p3(T0 


,0), 


then  a  process  of  lens  formation  will  occur. 
Finally,  if 


(4.130)  kfg( 0)  <  kfTp  +  (H-  b)LPi<p3(Tp,0) 
then  a  process  of  frost  penetration  will  take  place. 

The  diagram  on  figure  4.7,  where  we  put  the  initial  temperature  A(0)  on  the  *-  axis  and 
the  temperature  *(0)  on  the  t/-axis,  exhibits  the  fact  that  for  any  pair  of  values 

(fc(0),«,(0))  chosen  on  the  quarter  of  plane  A(0)  >  0,  ,(0)  <  0,  the  kind  of  process  is 
discriminated. 
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kfg(0) 


formation 


frost  penetration 


fig.  4.7:  the  straight  line  5,  corresponds  to  kfg(0)  =  kjT(r  +  (H-b)Lpi<p3(T^0),  the 


curve  S0  to 


M( 0)  =  W*( 0))  +  (B-  -  b)Lp,v3(Tp(h( 0)),0){1  + 

The  5p,m  curve  8iven  by  the  equation  kfg(0)  =  k }T p  +  (H  -  b)Lpitp3(Tp, 0) 
corresponds  to  Tp  =  Tm  (see  eq.  (4.4.58)).  The  two  regions  that  form  the  region  of  lens 
formation  correspond  to  the  cases  1)  and  2)  introduced  in  the  proof  of  proposition  4.1. 


>0- 


We  are  going  now  to  give  an  example  of  transition  from  one  process  to  the  other. 
We  choose  non  decreasing  boundary  temperatures 

(4.131)  h(t)>  0,  g(t)>  0 


and  we  assume  that  for  t  ^  0  (4.130)  is  verified.  The  process  of  frost  penetration  goes  on 
up  to  the  time  r  (see  remark  4.3)  when  (4.126)  holds: 


=  W)  +  W')  - 

where  Tp  is  the  solution  of  (4.102). 

Define  b1  =  zs(r),  H1  =  zT(r).  Consistently,  the  value  b  in  the  definition  of  the  functions 
‘Pi  and  <p3  has  to  be  replaced  by  6a: 


-  PwKo  S  ’  ^>(«,<)  =  J~r  ( J  kfJ(n)dv  -  kuh(t)\ 

kuKt)~  J  kff(rj)dr,  '  '  0  ' 

0 

The  temperature  Tp(h{t))  does  not  depend  on  b,  as  we  see  from  (4.26),  (4.27). 

Thinking  of  t  =  r  as  the  starting  point  of  a  second  process,  we  notice  that  condition 
(4.129)  holds:  more  precisely,  it  is 


kfg(r)  =  kfTp(h{ r))  +  (Ht  -  bjLp^T. p(A(r)),r)[  1  +  p  L 


K2{Tp{h{T))\ 

kl}(Tp{h(T))f 


By  virtue  of  (4.64),  the  temperature  Ts(t)  solution  of  (4.21)  is  exactly  Tp(h(r)). 

Choosing  zs(t)  =  zs(t)  =  6j  and  zT(r)  =  Hv  we  solve  again  the  system  (Stmp)  +  (C)  +  (A),  in 

the  same  way  as  we  described  in  section  4.1  and  we  get  a  first  process  of  lens  formation. 
Notice  that 


(4.132)  =0)  WW  >0)  <>r 

Indeed,  Ts(r )  verifies  (4.67)  and  for  t  >  r  the  temperature  Ts(t)  is  strictly  increasing  by 

the  assumption  (4.131)  (cfr.  proposition  4.6),  while  Tp(h(t))  is  non  increasing  (see  remark 
4.2). 
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The  first  equation  of  (4.132)  guarantees  that  the  water  flux  qw  is  continuous  for  t  =  r,  as 
we  see  from  (4.5). 

Owing  to  proposition  4.6  and  to  (4.131),  the  process  of  lens  formation  either  stops  in  a 
finite  time  tf  or  tends  to  a  stationary  profile;  in  the  first  case,  for  t  >  tf  a  process  of 
melting  will  occur;  in  the  second  case,  the  solution  tends  to  the  asymptotic  values 
defined  by  (4.72),  (4.75). 

Let  us  examinate  now  the  possibility  of  getting  the  inverse  process,  that  is  from  lens 
formation  to  frost  penetration.  Obviously,  we  have  to  drop  the  assumption  (4.131)  and 
choose  rapid  freezing  profile  for  the  boundary  temperatures.  Actually,  it  is  necessary  for 
the  temperature  T5(i),  solution  of  (4.20),  to  go  under  Tp(h(t))  (cfr.  case  2)  in  proposition 
4.5  and  example  (4.97)). 

For  t>tj,  the  solution  describing  lens  formation  is  meaningless:  indeed,  condition  (4.9) 
is  violated,  since  it  would  be  (see  remark  4.2): 

(4.133)  <  0,  for  t  >  tf 


where  b  represents,  as  usual,  the  base  of  the  lens. 

Thus,  the  solution  of  (4.20)  can  not  be  accepted  for  t  >  tf.  Let  us  look  for  solutions 
describing  a  frost  penetration  process  for  t  >  if ,  taking  t  —  tf  as  the  starting  time  and 
imposing 

(4.134)  dpj*s(*)£)  =  o,  for  t  >  tf. 


The  updated  initial  conditions  axe  (see  (4.83)): 

zs(tf)  =  b,  zT{tf)  =  H  -  J  ^(Ts(r),r)dT 


In  the  definition  of  and  <p3  the  value  for  b  has  not  to  be  changed,  unlike  the  previous 
transition  example,  since  during  the  lens  formation  the  base  of  the  ice  layer  keeps  at 
rest. 

Calling 


(4.135)  (?L/)(t),  $\i),  pV\zs(i),t)) 
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the  solution  of  (4.1)-(4.5)  obtained  by  imposing  (4.134)  (in  particular,  =  Tp(h(t )), 

cfr.  (4.102)),  we  are  going  to  show  that  the  boundary  #>(*)  can  not  be  non-decreasing. 
Indeed,  let  us  denote  by 

(4-!36)  ($(<),  z<P(t),  z$(t),  T(!)(t),  p(J*(zs(t),t)) 

the  formal  solution  of  (4.1)-(4.5)  for  t  >  t}  obained  by  imposing  z^(t)  =  o  (thus  z^(t)  =  b ) 
that  is  not  a  solution  of  (Stmp)  +  (C)  +  (A),  because  (4.9)  is  violated  (see  (4.133)). 

Assume,  contrary  to  our  claim,  that  there  exists  an  interval  [tf,ia)  such  that 

(4.137)  4/}(<)>0,  te[tf,ta). 

We  have,  according  to  (4.133): 

(4.138)  T^(t)  <  Ty\t)  =  Tp(h(t)),  t  e  ( tf1ta) . 

Furthermore,  (4.12),  (4.99)  and  (4.137)  yeld 

(4.139)  4°(<)  =  b  <  zfiit)  >  zfl(t),  t  e  (tf,ta). 

From  (4.5),  (4.99)  and  (4.139)  we  get: 

(4-140)  ,«(<)>  te(tf,ta). 

Using  (4.3),  (4.137),  (4.140)  and  the  comparison  theorem  for  ordinary  differential 
equations,  one  has  (assume  pi  <  pj): 

(4.141)  4%)  <  te(tf,ta). 

But  (4.2),  (4.138),  (4.139)  and  (4.141)  entail: 

Qw\t)  >  9w\t),  te(tfjta). 

which  contradicts  (4.140).  We  deduce  that  #>(*)  can  not  be  non-decreasing. 
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Besides  that,  we  have  that  x(t)  =  y{t)  =  -  z^\t)  is  the  solution  of  system  (SF) 

defined  in  par.  4.2.  From  the  equations  of  that  system  we  easily  get: 

=('4«)fg+ m -m- £<<$§)*’<<). 

Since  (cfr.  (4.117)  and  (4.126)) 

we  can  state  that  the  solution  x(t)  can  not  oscillate  infinite  times  so  that  the  sign  of  x  is 
not  defined  in  the  neighbours  of  tj  t>tp  provided  that  the  given  temperatures  h  and  g 
are,  as  it  is  natural,  piece- wise  increasing  or  decreasing. 

According  to  the  fact  that  z^\t)  is  not  non-decreasing,  as  we  proved,  we  conclude  that 

4/}(<)  <  o 

in  a  suitable  interval  (<y,7). 

Hence,  condition  (4.98)  is  verified  and  the  solution  (4.135)  describes  a  frost  penetration 
process. 
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